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I.  INTRODUCTION 

The  study  of  an  important  and  challenging  problem  in  science  and  engineering 
has  been  the  understanding  of  the  strength  and  fracture  behavior  of  stressed  solid 
systems.  This  is  particularly  true  in  the  failure  behavior  and  its  prediction  in 
viscoelastic  material  systems.  Both  the  structural  and  functional  application  of  these 
materials  demand  a  better  understanding  of  their  behavior  and  failure  mechanisms. 
When  sufficiently  large  tensile  stresses  are  associated  with  these  materials,  various 
modes  of  failure  develop.  To  elucidate  these,  one  common  mode  of  response,  namely 
crazing  under  an  applied  simple  stress,  must  first  be  understood.  Major  advances  and 
breakthroughs  in  the  crazing  behavior  in  microscopic  and  macroscopic  levels  of 
understanding  will  yield  tremendously  useful  information  not  only  theoretically  but  also 
practically.  Considerable  technological  and  scientific  significance  is  attached  to  this 
proposed  endeavor.  The  initiation  and  propagation  of  crazing  as  quasifracture,  the 
time  dependent  fracture  strength  of  oriented  polymers,  the  associated  molecular 
orientation  and  ultimate  strength  in  and  around  a  craze,  the  interaction  of  crazes  in 
polymeric  and  composite  systems  are  just  some  of  the  features  to  be  understood.  The 
determination  of  the  time  dependent  fracture  strength  of  polymers  and  composite 
systems,  the  displacement  field  and  the  stress  distribution  in  the  vicinity  of  craze-crack 
transition  region  as  well  as  the  propagation  behavior  of  craze  and  crack  are  important 
problems  to  be  solved  prior  to  the  consideration  of  many  other  relevant  topics. 
Currently  a  firm  foundation  has  been  established.  It  appears  that  continued  research 
in  the  relevant  outgrowth  topics  will  result  in  a  truly  fruitful  understanding  of  the  subject 
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II.  BACKGROUND  INFORMATION  AND  OBJECTIVES,  with  REFERENCES 

Advanced  reinforced  plastics,  consisting  of  a  polymer  matrix  and  fibres, 
continue  to  generate  great  interest  in  their  application  to  high  performance  structural 
components.  Fracture  of  these  composite  systems  may  result  from  flaws  in  fibres  or 
matrix  as  well  as  the  failure  of  the  bonds.  Thus  the  strength  of  any  such  composite  is 
governed  by  the  time  dependent  strength  characteristics  of  the  matrix,  the  fibres  and 
the  bonds.  Because  of  the  difference  in  the  mechanical  behavior  of  the  three 
constituents  of  composites,  up  to  now  many  strength  criteria  have  been  considered 
and  developed  by  scientists  and  engineers  all  over  the  world  as  reflected  by,  for 
example,  several  recent  references  [1  through  3]. 

The  studies  of  time-dependent  failure  of  composites  have  been  relatively  scarce 
in  spite  of  the  strong  dependence  of  the  failure  characteristics  on  time.  The 
formulation  of  the  models  must  now  be  based  upon  the  microstructural  peculiarities  of 
deformation,  the  molecular  orientation,  temperature  and  time.  Aside  from  the 
phenomenological  models,  perhaps,  statistical  models  should  also  be  considered 
concurrently  so  that  they  may  reinforce  each  other's  findings  and  development. 

REFERENCES 

1 .  Handbook  of  Composites  Series  (Strong  Fibres,  Structure  and  Design,  Failure 
Mechanics  of  Composites  and  Fabrication  fo  Composites).  Edited  by  A.  Kelly  and  Yu. 
N.  Rabotnove,  North  Holland  (1985). 

2.  Proceedings  of  International  Symposium  on  Composite  Materials  and 
Stmctures.  Edited  by  T.  T.  Loo  and  C.  T.  Sun,  Beijing,  China  (June,  1985). 

3.  M.  F.  Kanninen  and  C.  H.  Popelar.  Advanced  Fractures  Mechanics.  Oxford 
University  Press,  New  York;  Clarendon  Press,  Oxford  (1985). 
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ill.  MICROMECHANICS  OF  POLYMERS  AND  COMPOSITES 

The  phenomenon  of  crazing  and  its  relation  to  some  fracture  analyses  are 
considered  as  follows: 

The  formation  of  a  craze  comes  about  from  a  physical  transformation  in  the 
deformation  processes  of  the  microscopic  material  molecules  under  tensile  stress. 
The  transformation  takes  place  from  a  homogeneous  deformation  to  a  craze 
configuration  when  a  critical  condition  is  reached.  Subsequently,  the  craze  boundary 
propagates  as  a  function  of  applied  stress,  time,  temperature,  physical  and  chemical 
influence  as  well  as  the  actual  microstructural  changes  subjected  to  geometrical 
constraints.  As  a  result,  usually  minute  voids  are  generated  among  oriented 
molecules  and  the  density  of  the  medium  in  the  crazed  region  is  nonuniformly  reduced 
whereas  the  bulk  of  the  homogeneous  material  body  deforms  more  uniformly.  The 
interface  boundary  layer  enveloping  crazes  of  many  solid  materials  is  capable  of 
being  drawn  and  transformed  into  bundles  of  highly  oriented  molecular  domain 
structure  in  the  craze  region.  Further  stressing  will  eventually  initiate  craze-crack 
transition.  It  appears  necessary  to  take  these  physical  variations  into  consideration  in 
any  mathematical  modeling  and  formulation  in  analyzing  the  stresses  from  the  time 
when  crazes  incept  to  the  time  when  they  propagate  and  transform  into  real  fractures. 

The  science  of  crazing,  a  quasifracture  state,  and  subsequent  cracking,  a 
fracture  state,  of  solid  material  systems  under  tension  has  been  making  large  strides  in 
the  recent  past.  The  crazing  mechanism  has  been  associated  with  molecular 
orientation  and  fracture  strength  [1  to  5).  Subsequently  various  methods  have  been 
utilized  to  determine  and  confirm  the  molecular  mechanism  with  respect  to  craze 
formation  and  fracture  in  thermoplastics  [5,6].  Essentially  under  tensile  stresses 
certain  solid  materials  deform  from  sites  where  high  stress  concentrations  are  created 
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and  crazes  develop.  Because  of  geometrical  constraints  and  energy  requirements, 
the  material  molecules  orient  themselves  in  the  direction  of  stressing  with  voids  among 
them.  As  stated  earlier,  the  presence  of  oriented  polymeric  molecules  in  a  craze 
region  bounded  by  surprisingly  smooth  interface  layers  is  visualized  as  an  actual 
physical  phase  transformation  in  the  deformation  processes  from  one  orientation  state 
to  another  depending  upon  the  magnitude  and  rate  of  applied  tensile  stress  [8], 
material  characteristics  as  well  as,  of  course,  temperature  and  physical  and  chemical 
environments,  etc.,  surround  the  solid  body.  As  a  result,  the  mechanical  behavior  of 
the  material  is  greatly  affected  by  the  macroscopic  geometry  and  the  distribution  and 
interaction  of  the  individual  crazes  as  well  as  the  microscopic  molecular  configuration 
and  voids  within  each  craze  region  and  along  its  immediate  boundaries  enveloping 
the  area.  Macroscopically  the  development  of  crazes  and  their  distribution  can  be 
detected  statistically  by  laser  diffraction  technique  [9].  The  geometry  of  an  individual 
craze  which  can  be  studied  by  focused  laser  beams  [10]  is  of  primary  importance  in 
understanding  the  processes  of  its  initiation  and  propagation  as  well  as  the 
deformation,  quasifracture-fracture  transition,  and  eventually  the  fracture  behavior  of 
the  medium.  Knowledge  of  craze  initiation  and  geometry  helps  in  determining  the 
craze  displacement  field,  the  stress  distribution  and  the  craze-crack  transition  and 
propagation  under  load  [1 1 ,12].  An  eventual  understanding  of  the  true  mechanism  of 
molecular  strength  and  fracture  behavior  of  a  simple  solid  matrix  and  a  complex 
composite  system  can  be  obtained  if  fundamental  microscopic  information  is  utilized  in 
macroscopic  analyses. 

In  a  craze  the  highly  strained  molecular  bundles  act  as  boundary  tractions  with 
great  strength,  any  governing  mathematical  formulation  must  include  this  feature  for 
any  adequate  analysis.  Crazes  of  different  forms  and  properties  have  occurred  in 


Crazing  in  Polymeric  and  Composite  Systems 
Annua)  Technical  Report,  April  30, 1988 


5 


polymeric  materials  [13,14]  and  other  solid  systems  including  even  single  crystals  [15]. 
An  analysis  is  highly  desirable  and  may  be  useful  for  studying  general  solid  systems. 

Both  long-  and  short-range  programs  may  be  considered.  It  appears  fruitful  that 
emphasis  be  placed  on  the  study  of  micromechanics  of  individual  craze-crack 
transition,  the  source  of  failure  under  various  internal  and  external  stresses  for  the 
matrix  and  the  composite  systems. 

The  nature  of  the  stresses  in  and  around  a  craze-crack  transition  region  is  the 
key  to  the  understanding  of  the  morphology  and  nucleation  as  well  as  the  propagation 
of  crazes  and  cracks.  The  first  attempt  in  calculating  the  state  of  macroscopic  tensile 
stress  field  in  the  direction  of  the  applied  load  as  a  function  of  craze  length  has  been 
based  upon  a  model  with  an  assumed  craze  boundary  displacement  as  a  crack 
opening  in  an  infinite  elastic  sheet  [16].  The  stresses  were  calculated  as  though  the 
craze  were  a  continuum  and  the  craze  boundary  developed  no  stress  perpendicular  to 
the  direction  of  applied  stress.  The  solution  of  the  two-dimensional-homogeneous 
biharmonic  equation  for  a  semi-infinite  elastic  medium  due  to  the  application  of  an 
external  pressure  to  the  surface  has  been  used  [17,18].  This  implies  that  the  craze 
behavior  is  independent  of  the  craze  medium  [18,19]  under  stress.  The  solutions  were 
obtained  using  a  Fourier  transform  technique  [20]  or  a  complex  variable  method  of 
analysis  [19,21,22].  With  proper  assumed  boundary  conditions  the  latter  method  of 
approach  gives  probable  stress  and  displacement  fields  surrounding  a  craze.  A 
model  for  craze  growth  has  also  been  considered  with  the  creep  of  craze  material  as 
the  cause  of  craze  propagation.  The  craze  growth  was  found  to  be  linear  with  respect 
to  the  log  of  time  [19]. 

The  aforementioned  stress  analyses  have  been  made  essentially  on  the  basis 
of  the  classical  elasticity  theory  for  a  homogeneous  elastic  medium  with  either  an 
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assumed  stress  distribution  for  certain  portions  of  a  crack  without  considering  any  time 
dependency. 

The  development  of  crazing  is  not  only  a  function  of  stress  but  also  a  function  of 
time  [23,24].  Using  the  current  theory  and  by  taking  into  consideration  the  isotropic 
and  anisotropic  material  constants  the  mathematical  model  describing  the  crazing 
mechanism  have  been  successful  [25]. 
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V.  OUTLOOK  FOR  THE  SECOND  YEAR 
(Work  to  be  considered  for  publication) 

1 .  "Crack-Induced-Craze" 

The  analysis  of  the  crack-induced-craze  in  polymers  is  believed  to  be  a 
fairly  general  phenomenon  in  fracture  studies.  This  work  will  deal  with  the  use 
of  a  viscoelastic  boundary  element  method  for  analyzing  a  polymer  quasi¬ 
fracture.  A  time  dependent  boundary  stiffness  will  be  considered  and  the 
viscoelastic  solution  in  the  time  domain  may  be  obtained  by  applying  the 
collocation  Laplace  inversion  technique.  Using  these  methods,  the 
quasifracture  problem  with  time  dependent  stiffness  fractions  in  a  two- 
dimensional  case  may  be  analyzed.  Both  the  craze  opening  displacement 
profile  and  the  envelope  stress  distribution  around  a  craze  can  be  computed. 
This  will  pave  the  way  in  evaluating  the  propagation  history  of  both  the  crack 
and  the  craze.  Results  thus  obtained  may  be  compared  with  those  obtained  by 
previous  considerations  such  as  the  use  of  Dugdale  model  and  the  concern  on 
the  stress  concentration  phenomenon. 

2.  "Crazing  as  Quasifracture" 

Before  any  real  fracture  develops  under  stress  in  polymeric  or  composite 
systems,  it  seems  that  a  fairly  general  picture  common  to  most  solid  systems, 
crazing  incepts  first.  Following  the  previously  stated  craze-crack  transition  and 
crack-induced-craze,  the  initiation  of  crazing  is  simply  a  special  case.  In  the 
absence  of  crack  the  craze  as  quasifracture  has  been  studied  by  many 
scientists.  Since  the  boundary  element  method  has  become  recently  a 
powerful  technique  for  solving  boundary  value  problems  including  some 
nonlinear  ones,  it  is  especially  important  as  a  tool  to  be  used  in  problems 
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having  viscoelastic  deformations  and  fractures.  Therefore,  it  may  be  fruitful  in 
developing  proper  procedures  for  calculating  the  stress  distributions  around  a 
craze  envelope. 

3.  Time  Dependent  Fracture  Strength  of  Solid  Bodies" 

Statistical  theories  in  fracture  kinetics  constitute  a  very  important  role  in 
investigating  the  fracture  strength  of  solids  and  their  utilization  in  modern 
engineering.  In  this  short  report,  a  review  of  some  of  the  recent  concepts  and 
models  is  provided.  The  main  concern  is  the  effect  of  the  breaking  stress  on  the 
time-to-break.  Based  upon  the  consideration  of  the  fraction  of  integrity  of  a 
medium,  a  number  of  models  have  been  evaluated  and  compared.  Two  basic 
considerations  used  for  evaluation  and  comparison  are  Zhurkov's  empirical 
kinetic  relationship  and  Hsiao's  statistical  absolute  reaction  rate  model.  Other 
considerations  reducible  from  these  two  are  also  given  for  comparison.  Using  a 
well-known  numerical  analysis  method,  it  appears  that  the  nonlinear 
mathematical  consideration  is  more  realistic  in  describing  the  time-dependent 
fracture  strength  behavior  of  a  medium  over  any  linear  ones.  The  computed 
results  seem  to  fit  reasonably  well  with  the  general  observations. 

VI.  OUTLOOK  FOR  THE  THIRD  YEAR 
(Research  to  be  conducted) 

1.  Temperature  Variation  during  Polymer  Failure" 

This  paper  attempts  to  discuss  the  temperature  variation  during  polymer 
failure  using  a  statistical  absolute  reaction  rate  theory.  At  fracture,  the 
temperature  may  increase  or  decrease  depending  upon  a  quantity  named 
fraction  or  integrity  /  and  its  rate  /  and  accelerator  /  as  well  as  a  stress  modifier 
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For  over  a  century,  scientists  and  engineers  have  observed  temperature 
variations  during  loading  and  testing  of  solids.  Most  work  in  this  area  focused 
on  metallic  systems:  temperature  changes  during  elastic  and/or  plastic 
deformations,  as  well  as  theoretical  investigations  based  on  mechanics  and 
thermodynamics. 

Using  the  statistical  absolute  reaction  rate  theory,  the  present  work 
attempts  to  analyze  the  temperature  variation  during  polymer  failure. 

2.  "Irradiation  on  Crazing  Development" 

The  aerospace  industry  of  the  mid  '80s  is  now  being  heavily  involved  in  an 
unprecedented  number  of  space  projects  which  require  the  use  of  a  diverse 
range  of  devices,  all  of  which  are  exposed  to  much  higher  levels  of  radiation 
than  the  ones  inside  the  earth's  atmosphere.  Polymeric  materials  are  in 
common  use  although  it  is  known  that  changes  in  their  molecular  structure  are 
induced  when  exposed  to  high  energy  irradiation.  The  effects  are  not  the  same 
for  all  polymeric  systems  since  some  of  them  may  degrade,  others  are  induced 
to  crosslinking  and  still  others  may  crosslink  to  a  saturation  point  and  then 
degrade.  In  all  cases,  the  mechanical  properties  are  affected  which  implies 
changes  on  the  performance  and  longevity  of  the  material  system. 

At  present  there  is  a  need  to  understand  the  mechanism  responsible  for 
these  alterations  as  well  as  a  model  from  which  analysis  and  prediction  of  the 
polymer  crazing  behavior  in  environments  exposed  to  radiation. 

3.  "Three-Dimensional  Crazing" 

Polymers  and  polymeric  composites  usually  fail  by  first  developing  crazing 
on  the  surface  of  the  material  system.  Internal  crazes  can  also  be  initiated  when 
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sufficient  and  necessary  conditions  exist.  The  time-dependent  craze  failure 
process,  whether  two-dimensional  or  three-dimensional,  may  be  characterized 
by  several  stages:  deformation,  development  of  microporosity,  craze  initiation, 
craze-crack  transition  and  propagation  until  complete  failure  occurs.  The 
interrelationship  among  the  applied  stress,  craze  initiation,  time  and 
temperature  has  been  established  and  a  fairly  general  time-dependent  theory 
on  craze  initiation  in  viscoelastic  media  has  been  formulated.  This  craze 
initiation  criterion  is  three-dimensional.  It  would  be  highly  desirable  to  examine 
this  new  criterion  and  analyze  the  craze  in  three-dimensions. 


VII.  POSSIBLE  FUTURE  TOPICS  AND  THEIR  IMPACT  ON  POLYMERIC 
AND  COMPOSITE  RESEARCH 

"Analysis  of  Craze  Interactions" 

"Potential  Energy  Mimimization  of  Composite  Matrix" 

"Potential  Energy  of  Fiber  Reinforced  Composite  System" 

"Elasticity  Constants  of  Composites" 

"Thermal  Coefficients  of  Composites" 

The  aforementioned  individual  descriptions  are  some  immediate  projects  to  be 
completed  and  some  others  to  be  explored.  Whenever  possible  time  and  temperature 
dependent  stress  analyses  need  to  be  introduced  in  the  constitutive  modeling  of  the 
mechanical  behavior  of  composite  systems  at  micro-  and  macro-structural  levels.  This 
interest  is  in  the  development  of  analytical  methods  as  tools.  To  be  followed  are  some 
experimental  verifications  and  certain  temperature  data  for  quantitative  analyses  and 
numerical  computations.  It  is  hoped  that  these  projects  will  have  an  impact  leading  to 
a  better  understanding  of  the  effect  of  time  and  temperature  on  craze  and  fracture  of 
composite  systems. 
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VIII.  SUMMARY 

*  Noncontinuum  Craze-Crack  transition  yields  damage  behavior  of  composite 
systems. 

*  Studies  of  time  dependent  failure  of  composite  systems  are  important 
practically  and  theoretically. 

*  Objectives  are  to  obtain  better  understanding  of  the  connection  between 
microstructure  and  mechanics. 

*  An  overall  view  is  given  during  the  first  year,  in  the  published  paper  "Non- 
Continuum  Craze-Crack  Transition." 

*  During  the  second  year,  several  topical  papers  will  be  published.  They  are 
related  to: 

"Crack-Induced-Craze" 

"Crazing  as  Quasifracture" 

"Time  Dependent  Fracture  Strength  of  Solid  Bodies" 

*  Some  research  work  will  be  considered  for  the  third  year": 

"Irradiation  on  Crazing  Development" 

"Three-Dimensional  Crazing" 

"Temperature  Variation  During  Polymer  Failure" 

*  Possible  future  research  topics: 

"Analysis  of  Craze  Interactions" 

"Minimization  of  Potential  Energy  for  Composite  Matrix" 

"Minimization  of  Potential  Energy  for  Composite  Systems" 

"Analyses  of  Elasticity  Constants  of  Composites" 

"Analyses  of  Thermal  Coefficients  of  Composites" 

*  Six  (6)  copies  of  "Non-Continuum  Craze-Crack  Transition"  are  appended. 


NON-CONTINUUM  CRAZE-CRACK  TRANSITION 
C.C.  Hiiao 
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ABSTRACT 

An  approach  studying  the  damage  problem  ot 
polymeric  and  composite  material  systems  is  reviewed,  it 
appears  highly  desirable  to  investigate  the  strength  and 
fracture  behavior  of  stressed  sold  systems  by  combining 
microstructural  information  with  macro-  mechanical 
analysis. 

Starting  from  the  initiation  of  crazes,  the  process  of 
transition  of  a  craze  to  crack  in  an  infinite  viscoelastic 
medium  under  stress  is  described.  The  actual  physical 
change  of  a  craze  into  a  crack  is  considered.  Emphasis  is 
placed  on  the  nature  of  its  time  dependency.  The 
enlargement  of  both  craze  and  crack  is  analyzed  under  a 
simple  state  of  tension.  Suggestions  on  future  research 
upon  temperature  influence  on  crazing  and  craze-crack 
interactions  are  also  given. 

NOMENCLATURE 

A.B  material  constants 

an  integer  constants  (-<»  <  n  <  °°) 

a(t)  time-dependent  length  to  high  stress  level  cr2 

b(t)  time-dependent  crack  length  measured  from 

center  of  symmetry 

c(t)  time-dependent  craze-crack  length  measured 

from  center  of  symmetry 

Cb(t)  time-dependent  material  creep  compliance 

function 

6f  (t)  rate  of  energy  dissipation  of  creue  fibril 
domains 

d(x,t)  time-  and  position-dependent  diameter  of 

craze  fibril  domains 


d  t  average  diameter  of  fibril  domain  under  craze 

envelope  stress  O  \ 

d2  average  diameter  of  fibril  domain  under  craze 

envelope  stress  o2 

E(t)  relaxation  modulus 

ejj  deviatoric  strains 

f  fraction  of  integrity  of  microstructural  system 

Gi(t-r)  time-dependent  deviatoric  relaxation  modulus 

G2(t-r)  time-dependent  dilatational  relaxation 
modulus 

Hc(t)  rate  of  energy  absorption  of  quadrantal  craze 
h,l2  constants 

J2(£-t\  )  time  dependent  bulk  creep  compliance 
function 

K),  K2,  K3,  Kg  differential  or  integral  functions 
Kb  -  o)bexp(-  U/RT  *fif( t)J  rate  coefficient  of  breakage 
Kr  ■  torexpj-  U/RT  -  V'f'(t)]  rate  coefficient  of  reformation 
Kf(t)  rate  of  kinetic  energy  of  craze  fibril  domains 
n  integers 


nj.nj  unit  vectors  (i.  |  *1,2  or  3) 

P(t)  time-dependent  load 
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dP 


infinitesimal  elemental  vector  in  reference 
frame  X«  (K  - 1 , 2  or  3) 


z  =  x3 


third  coordinate  axis 


«0 

constant  linear  thermal  coefficient  of 

dp 

infinitesimal  elemental  vector  in  current  frame 

Xh  (k-1.20f  3) 

* 

expansion 

positive  parameters 

P 

craze-crack  system  depth  measured  from 

center  of  symmetry 

h 

delta  function 

R 

universal  gas  constant 

n>,  rc 

material  constants 

S 

magnitude  of  deviatoric  stress  tensor 

£ 

orientation  strain  {- 1  <  £  <  ««) 

S 

deviatoric  stress  tensor 

<kl 

strain  tensor  (k.  1-1.2  or  3) 

s 

Laplace  parameter 

6.4 

spherical  coordinates 

Sii 

deviatoric  stresses 

0(T) 

temperature  function 

'1 

S(t) 

rate  of  creation  of  craze  fibril  domain  surface 

X(x,t) 

time-  and  position-dependent  draw  ratio 

T 

absolute  temperature 

Xt 

draw  ratio  beyond  craze  mid-section 

t 

real  time 

under  stress  d , 

fa 

time  when  a(t)  reaches  position  x 

Xml 

draw  ratio  of  craze  mid-section  at 
stress  level  d , 

tb 

time-to-break 

Xm2 

draw  ratio  of  craze  mid-section  at 

*c 

craze  initiation  time 

stress  level  d2 

tf 

craze-crack  transition  time 

V 

strain  ratio 

tn 

time  when  b(t)  reaches  position  x 

S-t+fT), 

,  t\  -  r4(T)  shift  times 

fn 

nucleus  incubation  time 

n 

strain  ratio 

tx 

time  for  tip  of  craze-crack  system  to  reach 

P 

density  of  the  probability  distribution 

position  x 

function  of  molecular  orientation 

h 

period  of  time  fibril  domains  subjected  to 
stress  d2 

a 

appied  stress 

O0(T.t) 

constant  temperature,  time-dependent  stress 

u 

Uf(t) 

activation  energy 

<*I.<J|| 

principal  stresses 

rate  of  strain  energy  absorption  of  craze  fibril 
domains 

cf  i,  d2 

craze  envelope  stress  levels 

Vf 

volume  fraction  of  craze  fibril  domains 

db 

tensile  strength  or  fracture  strength 

u',i 

displacement  gradients 

dc 

craze  envelope  stress 

w(x,t) 

craze  opening  displacement  measured  from 

dij 

stress  tensor  (i,i>i.2or3) 

Wo(*.t) 

center  of  symmetry 

*uz(xAt).  |x|<c(t).  opening  displacement 

Oii 

time  derivative  of  the  isotropic  stress  tensor 

measured  from  center  of  symmetry  of 
craze-crack  system  as  defined 

<*max 

maximum  breaking  stress 

X  *  X, 

coordinate  along  first  coordinate  axis 

Omin 

minimum  breaking  stress 

XK 

coordinates  In  reference  frame  (K»  1,2  or  3) 

r 

time  or  dummy  parameter 

Xfc 

coordfriates  in  currant  frame  (k  ■  i,  z  or3) 

♦(x.t) 

stress  function  at  spatial  position  x  and  time  t 

$(T)  temperature-tune  shift  function 

'f  axial  stress  of  an  efement 

d<o  soM  angle 

6>b  frequency  of  motion  in  breaking  process 
u>r  frequency  of  motion  in  reforming  process 


PRELIMINARIES  AND  INTRODUCTION 

For  a  long  time,  the  analyses  and  predtetion  of  the 
mechanical  strength  behavior  of  engineering  components 
have  been  dependent  mostly  upon  the  application  of 
continuum  theories.  This  is  true  even  in  failure  studies. 
Elasticity,  viscoelasticity,  plasticity  and  linear  elastic 
fracture  mechanics  have  dominated  much  of  the  analytical 
investigation  in  solid  mechanics.  Failure  criteria  for 
continuous  media  are  considered  to  be  independent  with 
respect  to  the  integrity  of  the  media.  There  is  no 
introduction  of  noncontinuum  damage  mechanics  based 
upon  microstructural  behavior,  and  discrete  interactions 
among  discontinuous  flaws  are  not  predictable. 

Recently,  however,  a  new  research  direction  has 
emerged  to  connect  microscopic  material  behavior  with 
macroscopic  structural  mechanics.  This  is  motivated 
primarily  by  the  strong  desire  to  design  the  constitution 
and  configuration  of  the  microstructure  of  material  systems 
for  obtaining  required  macromechanical  properties  and 
functions.  In  doing  so,  the  integrity  of  the  media  is 
included  in  failure  criteria  as  well  as  in  the  constitutive 
description.  Micromechanisms  and  their  interactions  are 
taken  into  consideration  in  analyses  aside  from  their  time 
and  temperature  dependencies. 

In  failure  investigations  of  solids  and  composite 
material  systems,  an  Important  and  challenging  problem 
in  science  and  engineering  has  been  the  attempt  to 
understand  the  strength  and  fracture  behavior  of  stressed 
solid  systems.  Either  continuum  or  noncontinuum  method 
of  approach  has  been  used  to  study  the  problem.  The 
continuum  damage  mechanics  approach  deals  with  the 
phenomenological  behavior  of  matrix  cracking  in  fiber- 
reinforced  composites  using  quantities  such  as  stress, 
strain,  strain  rate  and  temperature  fields.  The 
noncontinuum  approach  deals  with  the  individual  micro¬ 
damages  such  as  minute  crazes  and  cracks  formed  in  the 
matrix  together  with  the  matrix-fiber  interfaces  in 
composites.  Usual  field  quantities  appropriate  to  the 
problem  as  wen  as  unusual  parameters  are  introduced  to 
effect  a  satisfactory  solution.  During  such  a  course  of 
investigation,  both  micro-  and  macro-structural  information 
is  connected.  This  noncontinuum  approach  reties  on  the 
mode  of  damage  and  micromechanism.  Aside  from  usual 
parameters,  the  ultimate  results  can  depend  upon  other 
quantities  such  as  a  function  of  the  fraction  of  integrity 
andtor  a  distribution  function  of  molecular  orientation,  ole. 
The  following  gives  an  example  concerning  a 
noncontinuum  micromechanics  and  craze-crack  transition 
behavior. 

Since  the  first  pubtication  of  the  pioneer  work  on 


macroscopic  behavior.  However,  little  attention  has  been 
given  to  make  this  connection  for  many  decades. 

One  of  the  first  analyses  of  the  growth  of  crazing 
was  modeled  as  a  continuum  theory.  Under  a  critical 
tensile  stress  the  linearly  elastic  material  initiates  crazes 
which  are  rate  insensitive.2  Several  other  analyses  using 
elasticity  theory  have  also  appeared  more  recently.3-4 
Similar  to  the  analysis  of  a  craze,  many  craze-related 
crack  problems  have  been  investigated  in  viscoelastic 
media.5'10  The  intention  of  these  studies  was  to  take 
care  of  the  energy  dissipation  which  takes  place  in  the 
viscoelastic  bodies  not  considered  in  the  elastic  theory. 
The  time-dependent  size  and  shape  of  cracks  in  linearly 
viscoelastic  isotropic  continuum  media  have  been 
analyzed.  No  microstructural  nature  is  considered 
however.  Subsequently,  approximate  methods  of 
analysis  have  been  presented  and  illustrated  with  a 
failure  zone  to  obtain  viscoelastic  stresses  and 
displacements  for  elastic  solutions.  In  the  analysis,  in 
order  to  satisfy  the  developed  governing  equations,  the 
failure  zone  in  the  elastic  problem  is  modified  to  meet  the 
traction  boundary  condition  for  the  crack  faces.  In 
general,  the  connection  between  the  mechanics  of  the 
problem  and  the  noncontinuum  microstructural  nature  is 
not  made  in  these  investigations.  Therefore,  in  the 
following  example  a  noncontinuum  craze-crack  transition 
is  analyzed  in  light  of  the  actual  microstructural  behavior 
and  the  mechanics  of  the  problem.  Equations  governing 
the  rupture  of  the  fibril  domain  structure  in  the  middle  of 
the  craze  envelope  surfaces  subject  to  a  uniform  tension 
is  considered  for  an  isolated  craze  in  an  infinite « 
viscoelastic  polymer  sheet.  Solutions  yield  both 
information  on  the  time  dependent  craze-crack  transition 
and  the  displacement  field  around  the  craze-crack 
envelope  profile  describing  the  shape  of  the  craze-crack 
region.  Before  this  is  done  relevant  noncontinuum 
information  is  reviewed.  Justifications  and  significance  of 
using  the  noncontinuum  microstructure  are  described. 


TIME  DEPENDENT  DEFORMATION  AND 
MOLECULAR  ORIENTATION 

By  incorporating  microstructural  information,  the 
deformation  of  a  material  system  may  be  analyzed  under 
stresses  with  the  help  of  classical  continuum  theory.  For 
some  polymeric  and  composite  systems,  the  deformation 
processes  may  be  characterized  to  contain  a  molecular 
orientation  mechanism  with  a  microporosity  sensitive  to 
time  and  temperature. 

A  realistic  medium  may  be  represented  by  a  system 
of  microstructural  elements  which  translate  and  rotate 
under  stressing.  As  shown  in  Figure  1 ,  an  elemental 
vector  dP  in  the  reference  frame  Xk  (K  »  1 ,  2,  or  3) 
transforms  to  dp  in  the  current  frame  x*  (k  « i.  2  or  3)  under 
a  time-dependent  load  P(t).  Then  the  stress  tensor  oy  at 
any  point  can  be  calculated  under  certain  conditions.11 


<Jij(«w,T.t)  •  Jp(e, feni)  He.fe,T,t) 


*(e,f€,T,t)nsnjdtt, 


crazing,1  there  was  an  inadvertent  attempt  to  connect  the 
microstructure  of  the  polymer  medium  with  Its 
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Fig.  1  Molecular  Entanglement  and  Orientation  m  Deformation 

where  <kl  is  the  strain  tensor, 

T  is  the  absolute  temperature, 

t  is  time, 

p  is  the  density  of  the  probability  distribution 

function  of  molecular  orientation, 

d.  <J>  are  spherical  coordtoates, 

f  is  the  fraction  of  integrity  of  the 
microstructional  system, 

f  is  the  orientation  strain  (- 1  <  t  <  «>), 
is  the  axial  stress  in  an  element, 

nj,  nj  are  unit  vectors,  and 

do  equals  sinOdOd^  representing  the  solid  angle 
from  which  statistical  expectation  may  be 
evaluated. 


Since  the  stresses  are  functions  of  the  orientation 
strain,  it  is  likely  that  the  constitutive  behavior  will  be 
greatly  affected.  A  simple  illustration  is  shown  in  Figure  2 


where  the  nonlinear  deformational  behavior  of  an 
isotropic  continuum  medium  must  be  modified  by  the 
molecular  orientation  behavior  to  give  a  more  realistic 
combined  constitutive  description.  This  no  doubt  will 
affect  the  further  response  of  the  material  system  to 
stressing  such  as  the  initiation  of  crazes  and  the  stability 
or  instability  of  the  system. 

ENERGY  STATE,  TEMPERATURE  ANO 
CRAZE  INITIATION 

At  a  constant  temperature,  the  stability  of  the 
homogeneous  deformation  of  a  real  material  system 
under  a  simple  uniform  tension  will  eventually  be  upset  at 
some  time  when  a  specific  energy  state  of  the 
microstructural  system  develops.  The  possible  responses 
may  be  described  with  the  help  of  Figure  3.  Under  load. 


TIME 


Fig.  3  Lowest  Free  Energy  State  and  Deformation  Behavior 

the  stability  or  instability  of  a  material  system  may  best  be 
determined  by  its  free  energy  state.  The  material  system 
readjusts  its  microstructural  configuration  to  maintain  its 
lowest  energy  state  until  certain  dominating 
characteristics  develop.  After  an  incubation  period  is 
reached  the  material  system  may  deform  in  shear 
primarily  as  a  result  of  the  rotational  motion  of  certain 
microstructures.  Thus,  a  simple  shearing  band  develops 
in  the  system.  The  material  system  may  develop  cracks 
as  a  result  of  a  basically  translational  motion  of  the 
microstructure  including  possibly  the  slipping  and 
rupturing  of  molecular  bonds.  This  may  be  the  simple 
cracking  of  a  material  system.  When  a  complex  kinetic 
situation  occurs,  both  the  rotational  and  translational 
motion  of  the  microstructure  may  take  place.  Crazing  may 
develop  as  a  result. 

For  a  two-dimensional  craze,  the  isocnronous 
biaxial  locus  for  craze  initiation  is  given  in  Figure  4. 
Various  criteria  put  forward  to  date  are  plotted  for 
comparison.  Detailed  information  can  be  found  in  an 
earlier  reference.12  Only  the  three-dimensional  craze 
initiation  criterion  is  given  below  in  a  series  form  with  an 
as  an  integer  constant  and  n  an  integer. 
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Fig.  S  Stress,  Temperature  Effects  on 
Craze  Initiation  in  Polymers 


Stress  in  MPa 

Biaxial  locus  for  criteria  of  craze  initiation  in  polymers  by: 

1.  Sternstein  and  Ongchin, 

2.  Bowden  and  Oxborough, 

3.  New  Criterion, 

4.  Argon, 

5.  Distortion  strain  energy,  and 

6.  45°  reference  line 


Fig.  4  Biaxial  Locus  for  Criteria  of  Craze  Initiation  in  Polymers 


S(xk,T,t)2  X  an[3<x0e(T) 
n--«> 

t 

♦  JJ2^Tl)<3ii(xk.t)dr]n  (2) 


where  S  is  the  magnitude  of  the  deviatoric  stress  tensor  S 
which  must  overcome  an  intrinsic  flow  resistance,  an  are 
constants,  is  the  constant  linear  thermal  coefficient  of 
expansion,  d(T)  is  the  temperature  function,  J2(£  -  tj) 
is  the  bulk  creep  compliance  function  with  £  -  tf  (T) 
and  i)  *  r$(T)  as  shift  times  defined  by  the  "temperature- 
time  shift*  principle  for  “thermo theologically  simple" 
viscoelastic  media  and  0u  is  the  time  derivative  of  the 
isotropic  stress  tensor.  This  three-dimensional  craze 
initiation  criterion  is  reducible  to  any  of  the  other  criteria  by 
introducing  appropriate  values  ter  an-  In  Figure  4  this  new 
criterion  is  represented  by  line  3  for  the  two-dimensional 
situation.  This  Hne  reduces  to  each  of  the  other  curves 
from  1  through  5  when  appropriate  values  of  and  n  are 
introduced. 

The  influence  of  temperature  motion  on  craze- 
initiation  is  also  plotted  in  Figure  5  together  with  the 
applied  stress  for  the  polystyrene  material  system 
computed  with  ^  ■  -  30  MNrtn2,  a-i  «  0.3  MNrtn*  and 

ai«a 


MICROSTRUCTURE  AND 
MECHANICAL  STRENGTH 

Based  upon  the  microstructural  behavior  and  the 
dynamic  nature  of  molecular  motion,  the  fraction  of 
integrity  can  be  determined  and  utilized  to  establish  a 
fracture  criterion  for  s  material  system.  Subsequently  the 
connection  between  the  microstructure  and  the 
macromechanical  strength  can  be  made.  This  has  been 
done  by  means  of  the  statistical  theory  of  the  absolute 
reaction  rate.  Not  only  can  the  microstructural  orientation 
and  the  rupture  of  the  microscopic  structural  units  be 
incorporated  but  their  reformation  can  also  be  included  in 
the  theory  to  obtain  the  time-dependent  mechanical 
strength  particularly  for  long  times  as  illustrated  in  Figure 
5.11,13-15  This  tensile  strength  which  is  the  fracture 
strength  as  a  function  of  time  has  been  well  established 
as  shown.  The  tensile  strength  of  a  stressed  solid 
material  system  for  both  short  or  long  times  tends  to  level 
off  as  shown.  This  means  that  the  tensile  strength 
becomes  independent  of  time  for  very  short  times  as  well 
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Fig.  6  Micrestructural  Orientation  and  Reformation  Effects  on  a 
Time -Dependent  Mechanical  Strength 


as  for  very  long  times.  In  between  these  times,  the  tensne 
strength  is  Hnearty  proportional  to  the  logarithm  of  time. 
The  temperature  effect  on  the  time-dependent  mechanical 
strength  is  given  in  Figure  7  to  show  that  temperature  is 
an  important  entity  in  this  complex  situation  which  must  be 
taken  into  consideration  in  the  analysis  for  any  material 
system. 


LOGARITHM  OF  TIME-TO-BREAK  In  ^ 


Fig.  7  Temperature  Effects  on  Time-Dependent 
Mechanical  Strength 


Consider  a  material  system  composed  of  a  large 
number  of  randomly  oriented  similar  molecular  units  or 
elements  which  elongate  and  rotate  in  a  stress  field.  In 
order  to  formulate  a  temperature-  and  time-dependent 
failure  criterion  for  such  a  material  system,  one  can  use 
the  quantity  f(t).  the  time-dependent  fraction  of  integrity. 
By  calculating  its  time  rate  of  change  as  follows: 

f-Kr(1-f)-Kbf  (3) 

where  Kr  -  <i>r*xp(-  U/RT  -  Jf'l'(t))  is  the  rate  coefficient  of 
reformation  of  the  disconnected  units  and  K&  - 
o>bexp{-  U/RT  ♦  £¥(!)],  the  rate  coefficient  of  breakage  of 
the  connected  elements,  <i>r  and  o>b  are  respectively  the 
frequencies  of  motion  with  respect  to  forming  and 
breaking  processes  of  these  units,  U  is  the  activation 
energy,  R  is  the  universal  gas  constant,  and  7  and  0  are 
positive  parameters  which  modify  the  true  axial  stress  H'ft) 
in  each  elemental  unit. 

The  failing  of  a  material  system  is  when  f 
approaches  zero.  In  these  formulations,  Kr  and  Kb,  are 
both  functions  of  temperature  and  the  true  stress  in 
individual  elements.  Once  a  stress  0  is  applied  to  the 
material  system,  the  energy  state  is  altered  and  the  time- 
to-break  tjj  can  be  calculated. 

The  fracture  strength,  the  statistical  mean  strength 
in  the  vicinity  of  any  point  in  the  system,  has  been 
analyzed  and  found  to  be  proportional  to  the  modulus  of 
indhridual  elements,  their  length  and  the  number  of  the 
elements  per  unit  volume.  For  an  oriented  system,  the 
fracture  strength  to  a  function  of  deformation  as  briefly 
reviewed  in  the  previous  section.  For  a  fully  oriented 
system11  under  a  constant  applied  stress  0,  we  may 
write 


*(t)f(t)«0.  (4) 

The  functional  relationship  between  ty  and  t  is 

'I'2  a  o 
♦  --“[KrO'^-Kb^].  (5) 

Integration  results  in  the  time-to-break  to-  With  zero  initial 
time 

f  d» 

**“  J'*'[Kb+Kr<1-'f/0)]  ‘  (6) 

Usually  near  fracture,  Kr  becomes  unimportant,  and 
the  time-to-break  may  be  approximated  to 

J  •  <7> 

which  can  further  be  simplified  to 

tb-Aexp(-B0)  (8) 

with  A  and  B  as  material  constants. 

Now  let  us  consider  the  craze  problem  as  shown  in 
Figure  8.  Under  a  temperature  and  time  dependent  stress 
0o(T,t)  a  three-dimensional  craze  may  develop  from  the 
surface  of  a  material  body.  In  general  the  microstructure 
on  the  surface  of  the  material  body  is  composed  of  a 
network  of  highly  oriented  fibril  domains  drawn  out  of  the 
envelope  profile  and  separated  by  cavities.  This 
combined  structure  propagates  along  the  surface  of  the 
material  body  and  penetrates  into  the  body  as  indicated 
respectively  by  c(T.t)  and  q(T,t)  which  are  measured  from 
the  center  of  craze  at  the  origin  0.  The  stress  at  any  point 
in  the  material  is  designated  as  0(x1ax2,x3, T,t)  and  the 
craze  envelope  stress  by  0c(X|,x2,X3,T,t)  as  shown. 


x3  <To(T,t) 


Fig.  8  Microstructur®  of  a  Three-Dimensional 
Craze-Crack  System 


A  strong  interest  has  been  the  determination  of  the 
displacement  field  since  it  is  not  easily  measured.  This 
noncontinuum  feature  is  easily  seen  in  Figure  9,  in  which 
the  randomly  oriented  microstructure  is  being  drawn  into 
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exists  a  corresponding  curve  between  tp  and  the  applied 
constant  stress  0.  When  the  applied  stress  varies  with 
time,  one  can  use  the  linear  summation  damage  rule 
which  states  that  the  same  number  of  damages 
accumulate  during  identical  time  intervals.  Thus  the  time- 
to-break  is  determinable  by  the  equation:1* 

/s §»-’  <91 

By  adopting  this  microstructural  behavior,  the 
opening  displacement  w(x,t),  the  craze-crack  system 
length  c(t)  and  the  crack  length  b(t)  in  this  craze-crack 
transition  problem  can  be  calculated.  To  facilitate  the 
possibility  of  obtaining  numerical  solutions,  a  three-step 
envelope  stress  distribution,  as  shown  in  Figure  10,  is 
introduced: 


highly  oriented  fibril  domains.  A  schematic  diagram  of  a 
two-dimensional  craze  is  given  in  Figure  10  to  show  the 
pertinent  quantities.  For  simplicity  in  sequel  xj  is 
replaced  by  z  and  xt,  by  x. 


■o  o  <  x  <  b(t) , 
Oc(x.t)  »  -a  |  b(t)  <  x  <  a(t) , 
.02  a(t)<x<c(t). 


(10) 


Fig.  10  An  idealized  Two-Dimensional  Craze-Crack  System 
with  Assumed  Step  Envelope  Stress  Distribution 


ANALYSIS  OF  CRAZE-CRACK  TRANSITION 

Considering  this  two-dimensional  model,  an 
analysis  can  be  mads  to  obtain  the  time  dependent 
displacement  field  and  the  craze-crack  lengths  by 
incorporating  the  microstructural  behavior  into  the 
macromechanics  problem. 

Starting  from  the  highly  oriented  fibril  domain 
structure  in  the  craze  region,  the  required  time-to-break  tp 
is  calculable  based  upon  a  perfectly  oriented  molecular 
system  under  a  constant  stress  a.11'15  A  general  curve  tb 
vs.  O  is  given  in  Figure  6.  The  maximum  applied  stress  ijs 
indicated  by  <Tmax.  With  reformation  processes 
considered,  the  tensile  strength  of  the  material  system 
begins  to  deviate  from  the  dotted  curve.  There  exists  a 
minimum  strength  value  Cmm  tor  which  lb  goes  to  infinity 
asymptotically.14  Therefore  for  a  given  medium  there 


It  is  hoped  that  this  assumed  step  stress  function  will  yield 
good  approximations  in  both  the  displacement  field  and 
the  lengths  of  the  craze-crack  system.1 

With  the  three-step  envelope  stress  function,  the 
stress  sustained  by  the  fibril  domains  in  the  thin  mid¬ 
section  where  under  certain  conditions,  failure  has  been 
observed  to  occur  more  often,  is 


a(t)  <  x  <  c(t) , 
b(t)  <  x  <  a(t) , 


(11) 


where  Xm2<72  is  the  true  stress  bom  by  the  fibril  domains 
and  Xm2  is  the  draw  ratio  at  the  stress  region  02. 
Similarly  Xmt0y  and  Xmt  denote  the  corresponding 
quantities  under  region  Ov 

As  discussed  earlier,  under  large  stresses,  the 
time-to-break  tp  can  be  obtained  by  first  dropping  Kr  in  (3) 
to  get  (8),  then  introducing  (1 1 )  in  (9),  one  can  obtain: 

tp  m  Aexp(-  BXmi0i) 


-  t2(exp(B(Xm202-Xn»iai)J  -  1)  (12) 


where  t2  is  the  period  when  the  fibril  domains  experience 
the  higher  envelope  stress  a2. 

To  review  briefly  the  time-dependent  viscoelastic 
problem,  the  opening  displacement  measured  from  the 
center  of  symmetry  of  the  craze-crack  system  in  a 
viscoelastic  polymer  sheet  can  be  obtained  by  using  the 
well  known  correspondence  principle  in  linear 
viscoelasticity. 

The  field  equations  are: 


(13) 
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f#<X.Z,t)  +  Uj,i(X,Z.t)]  , 

(14) 

t 

Sg(t)-  J  Gt(t-r)det|(r) , 

-00 

(15) 

t 

<Jij(t)  -  /  G2(t-t)d<jj(T) , 

•00 

06) 

where  cry  is  the  stress  tensor,  fy,  the  strain  tensor,  ujj, 
the  displacement  gracfients,  G  \  and  G2  are  respectively 
the  deviatoric  and  dilatational  relaxation  moduli  of  the 
original  bulk  medium,  and  the  deviatoric  stresses  and 
strains  are  respectively 

Sirft)  -  cJjj(t)  , 

07) 

«ij(t)  -  f|j(t)  -jStjrkkW . 

(18) 

with  i,  j  and  k  *  x  or  z  denoting  dummy  variables  in  two- 
dimensional  problems.  The  boundary  conditions  are 
described  as: 


<*zx  -  o . 

|x|  <  c(t) 

(19) 

<?zz  (*.<M)  -  <*c  (x.t) 

,  |x|  <  c(t) 

(20) 

<7zz(x.z.t)  -  <J0(t) ,  - 

Oxx(*.z.t)  -  o . 

as  (x2+z2)  -->  «•  . 

(21) 

Oxz(x^.t)  -o , 

The  opening  displacement  w0(x,t)  is  defined  by 

w0(x,t)».u2(x,o.t).  M<c(t).  (22) 

To  solve  this  viscoelastic  problem.  The  Laplace 
transform  with  respect  to  time  is  applied  to  aW  of  the  field 
equations  and  the  boundary  conditions.  The  solution  to 
the  transformed  equation  can  be  found  by  using  the  well 
established  complex  variable  conformal  mapping  method 
or  the  complex  variable  stress  fonction  method  for  elastic 
medium. 

The  time  dependent  solution  of  the  problem  is 
obtained  by  Laplace  inversion.  This  solution  is  valid  only 
if  the  boundary  conditions  are  independent  of  time,  i.e., 
';t),  b(t)  and  oc(x,t)  remain  unchanged.  These 
restrictions  can  be  removed  by  using  a  sequence  of 
loading  and  unloading  steps,17  which  yields, 


t 

Wo(x.t)  -  Co(0)*(x,t)  ♦  j  etrft-t A(x^)dr .  (23) 

o 

where 

♦(x.t)  «  OeVc5^ - 1 J  Octf.t)  in(cWc*-ST  )d* 


x 

*K  ln<*  +  V*Z-5Z  )d£ 

0 

2  C 

♦;f<*e(i.t)lnfc«,  (24) 


Cb(t)  -  L-1  [(2(2G,+G*))  /  (s2G|(G|+2G2))]  (25) 

with  L~ 1  designating  the  Laplace  inversion  and  barred 
quantities  being  in  Laplace  domains. 

If  the  strain  ratio  vis  constant.  Cb(t)  reduces  to 

Cb(t)  -  L 1  [(2(1  -  8*V*))  /  S*E )  .  (26) 

where  E  is  the  relaxation  modulus  of  the  original 
unoriented  buk  polymer  medium. 


Taking  into  consideration  the  thickness  of  the 
primordial  layer  from  which  the  fibril  domain  structure  has 
been  pulled  out,  the  actual  opening  displacement  of  the 
system  becomes 

t 

w(x,t)-Cb(o)k(x1t)*  J  ib(t-r)*(x,r)dr 
o 


where  tx  is  the  time  when  the  tip  of  the  craze-crack  system 
first  reaches  the  point  x,  X  is  the  draw  ratio.  Since  the 
midsection  is  relatively  thin  in  a  craze,  the  X  function  can 


betaken  as 

rxi- 

b(t)<x<a(t) 

X(x,t). 

a(t)<x<c(t), 

(28) 

where  X|  is  the  draw  ratio  outside  the  mid-section  under 
the  envelope  stress  a,.  The  opening  displacement 
finally  is  obtained  as 


t 

♦  J  db(t-T)k(x,r)drJ ,  tor  a(t)  <  x  <  c(t) ,  (29) 
0 

-  » 

W-ftVWW)  +  J  CbO-rWx.rjdr] 

♦  {^-~y>»(x4a).  torb(t)<x<aW,  (30) 

t 

w(x,t)  ■  Ct>(o)*(x,t)  ♦  J  do(t-r)#(x,T)dr 
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♦  *  *r7)w(x,w  *  5T7wt*’ w  • 


loro<x<b(t)>  (31) 

where  ta  and  th  respectively  denote  the  times  when  a(t) 
and  b(t)  arrive  at  the  point  x. 

CRAZE  AND  CRACK  LENGTHS 

With  regard  to  the  length  of  the  craze-crack  system 
and  that  of  the  crack,  suppose  that  the  fibril  domain 
nucleation  rate  at  the  craze  tips  is  proportional  to  the 
growth  velocity  c(t)  of  the  system.  Then  the  energy  rate 
required  for  the  growth  of  craze  tip  is  r^ft)  and  that  for 
the  crack  tip  is  rb6(t),  where  rc  and  rb  are  material 
constants.  Based  upon  the  assumption  that  the  rate  of 
energy  required  for  drawing  the  fibrils  out  of  the  craze 
surface  envelope  is  proportional  to  the  rate  of  creation  of 
the  new  fibril  domain  surface  S(t),  using  the 
proportionality  constant  T,,  the  following  local  energy 
balance  equation  can  be  written  as: 

Hc(t)  -  rcc(t)  +  rbb(t)  ♦  r,$(t) 

+  Ur(t)  +  Df(t)  +  KfO)  (32) 

where  Hc(t)  is  the  energy  absorption  rate  of  the  quadrantal 
portion  of  the  craze  system,  Cl|(t)  is  the  strain  energy 
absorption  rate  of  fibril  domains,  C>r(t)  is  the  energy 
dissipation  rate  by  the  craze  fibrils,  and  Kf(t)  is  the  rate  of 
kinetic  energy  due  to  the  motion  of  the  craze  fibril 
domains.  Of,  Of  and  Kf  are  negligible  when  compared 
with  the  other  rate  quantities  in  the  quasi-static  conditions. 
With  this  simplification  and  the  terms  expressed  by 
elementary  parameters  defined  earlier,  the  following 
equations  can  be  established.  Since  the  stress  each  fibril 
domain  sustains  is  XOc  which  equals  CeA/|  if  Vf  is  the 
volume  fraction: 


Hc(t)-  J  <Jc(x.t)  — dx , 


c 


5(t )«  J 

o 


4 


Vr(x.t)  3w(x,t) 
d(x,t)  9t 


dx , 


(33) 


(34) 


where  d(x,t)  is  the  time  and  position  dependent  diameter 
of  craze  fibril  domains.  Now  using  these,  Equation  (32) 
becomes  an  implicit  nonlinear  deferential  equation  of  the 
craze-crack  system  length  c(t)  and  the  crack  length  b(t)  as 
follows: 


J  (Oc^^^^dx-rci^oli.  (35) 
0 

This  implicit  differential  equation  has  three 
unknown  quantities,  a(t),  b(t)  and  c(t)  to  be  determined.  In 
order  to  solve  this  equation,  some  subsidiary  equations 
are  necessary.  One  of  them  is  that17  the  craze  envelope 


stress  <?c(x.t)  must  balance  the  applied  external  load 
corresponding  to  an  applied  simple  tension  <J0  to  ensure 
that  the  stress  field  within  thn  uncrazed  bulk  material  is 
finite  everywhere  for  all  times  t.  In  mathematical  form  this 
means: 


J 


<3c(x.t) 

Vc2(t)-X2 


dx  =|cr0 . 


(36) 


When  the  three-step  envelope  stress  function  (10)  is 
substituted  in,  it  yields: 


a(t)-c(t)sin(£ 


(37) 


The  other  equation  is  from  the  consideration  of  the 
nature  of  the  failure  of  the  craze  material,  which  obviously 
provides  a  relation  between  the  craze-crack  system  and 
crack  lengths, 


bO)  =  c(t-tb).  (38) 

It  should  be  noted  that  the  time-to-break  tb  is  spatially 
dependent  and  it  is  evident  that 

b(t)-c[t-tb(b(t))].  (39) 


RESULTS  AND  DISCUSSION 

Generally  speaking,  c(t)  and  b(t)  can  be  obtained 
as  functions  of  time  by  solving  Equations  (35),  (37)  and 
(39).  But  it  is  still  rather  complicated  because  of  the 
unusual  form  of  Equation  (39). 

If  the  craze-crack  system  and  the  crack  propagate 
steadily  without  drastically  change  in  their  propagating 
speeds.  Equation  (39)  can  be  simplified  to  the  following 
form 


b(t)»c(t)-^c(t),  (40) 

where  tb  has  been  evaluated  and  displayed  as  Equation 
(12),  ie., 

tb  -  Aexp(-BXm10,) 


- 12{  explBtX^Oz  -  Xn^cr , )]  - 1} . 

Usually,  the  distance  c(t)  -  a(t)  is  relatively  small  since 
it  is  associated  with  the  region  of  stress  concentration. 
Thus  t2  can  be  expressed  as 


t2  ■  (c  -  a)  /  c . 


(41) 


Using  the  envelope  stress  profile  proposed  earlier. 
Equation  (24)  turns  out  to  be 


.  bVc2-x4  -  x-V  cs-b4 

*a,XlnbV^  +  xV^b7 
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Introducing  the  notations 


-  (02-0|)aln 


Vcz-x*  -  V c2-az 


-Ojbln 


Vc*-xa  -V c^-b2 . 
Vcs*x2  i/t^-b2 


♦ac  =  J  <Kx.t)d>c , 
a 


Since  the  product  of  the  average  fibril  domain  diameter 
and  the  envelope  stress  has  been  found  to  be  constant, 
that  is,  d  |<J  j  -  d 2a2  *  constant,  thus  let 

rd,,  b(t)  <  x  <  a(t) , 


I-O)  ,  0(1)  <X  <  aw  , 

•t). 

Ld2 ,  a(t)  <  X  <  c(t) , 


and  noting  that  V(  -  1/X  and  Equation  (32)  then,  by 
substitution.  Equation  (35)  becomes 


4rs ,  ,  dw 


^-dx.rcc  +  Tbb.  (44) 


^ba  =  J  ^x-t)dx' 

and  substituting  them  into  (44),  we  have 


K,a(t)  +  K2b(t)  ♦  Kjc(t)  +  K*  =  o , 


3$ar  d$hA 

K,  =  Cb(0)[  (l2-l,)*(a,t)  ♦  l2-gf-+  It-gJ1  1  • 

k2  -  cb(o){  I2^-  i,i  *  (b.t)  +  ^  » -  r„ . 

if  r  >~wi  .  .  9*ba  .  _ 

K3»cb(o)(i2-^-+i,  9c  )-rc. 


Since  0(c,  t)  -  o,  obviously, 


J 


K4  *  Cb(o)(l2®ac  +  h®ba) 


►gf  j*(*.t)dx] 


/  Cb(o)«(x,t)dx 


c  t 

-  /  J  £b(t-t)®(x,r)drdx } , 

a  0 


J  ^dx-  {Cb{o, [b*(b.t) - a$(a,t) 


hdt  J4><x-,)dxl 

b 


J  Cb(0)*(x,t)dx 


a  t 

f  f  Ct>(t-r)*(x,t)drdx } . 
b  0 


+ 12  J  Ct>(t-r)®ac<lr  + 1,  J  Cb(t-r)®badr  ,  (53) 


with  1 1  and  l2  being  constants,  and 

1541 

i2  -  (q2- '4rs  )  -~T-  (55) 

^2^m2  Xm2-1 

The  explicit  forms  of  those  quantities  in  expressions 
of  ICs  are  given  in  the  Appendix. 


The  calculations  should  be  divided  into  two  steps. 
First,  the  original  craze  propagates  during  the  absence  of 
a  crack.  This  can  be  calculated  by  simply  setting  b(t)  -  0 
in  Equations  (38)  and  (44),  which  degenerates  into  the 
case  discussed  earlier.17  After  certain  time  elapses,  the 
fibril  domains  first  produced  in  the  mid-section  of  the  craze 
breaks  down  and  crack  commences.  Second,  the  crack 
comes  Into  play  and  the  Equations  (38),  (41),  (42)  and 
(44)  must  be  used  to  calculate  simultaneously  the 
propagations  for  both  the  craze-crack  system  and  the 
crack. 

Now  to  illustrate  the  changes  of  a  craze-crack  system 
in  polystyrene,  a  Voigt  sold  is  taken  as  an  example.  The 
material  properties  are  taken  to  be1*-*3 
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J(t)  -  {29  +  5.0911  -  exp{-t)]  ♦  2.32[1  -  exp(-t/10)] 

♦  6.59  {1  -  exp(-t  /102)1  ♦  12.72J1  -  exp(-t  /10*)1 
+  0.71J1- exp{-t  /105)J  +  14.71  {1  -  exp(-t/106)J 


+ 1.02  [1-  exp(-t  /107)]}  •  10-s  m2/MN , 


A  -  exp(ll)  secs, 
B  -  0.05  m2/MN , 
X,  -  2.0, 

Xml  ■  2.5, 

Xm2  “  4 , 


d,  -20nm , 
d2  -  10  nm  , 

T,  -0.125  J/m2 , 
Tc  -  0.085  J/m2, 
r„-300  J/m2 . 


v  for  polystyrene  is  0.395,  a  constant,  and  the  applied 
stress  is  considered  to  be  36  MN  /  m2,  that  is:  O0  »  36 
MN  /  m2.  In  addition,  based  upon  some  experimental 
evidence,  <J,  and  o2  are  assumed  to  be  respectively 
40  MN  /  m2  and  60  MN  /  m2  for  numerical  calculations. 


complex  phenomena  are  yet  to  be  considered  in  the 
future  analyses. 


The  normalized  length  of  the  craze-crack  system 
and  that  of  the  crack  are  shown  in  Figure  1 1 .  where  the 


Fig.  1 1  Ome-Dependant  Crack  and  Craze-Crack 
System  Lengths 

normalization  is  made  with  respect  to  the  initial  length  of 
the  craze-crack  system.  The  opening  displacements  at 
different  times  are  plotted  in  Figure  12. 

For  ease  of  visualization,  the  normalized  craze 
length  as  a  function  of  time  is  also  given  in  Figure  13.  The 
non-linear  nature  of  the  craze  behavior  is  self-evident. 

As  can  be  seen  from  the  calculations,  both  the 
velocities  of  the  craze-crack  system  and  that  of  the  crack 
become  larger  and  larger.  Thus  further  calculation  is  not 
accurate  since  the  unsteady  propagation  gives  rise  to 
irregular  growth,  branching  and/or  bifurcation.  These 


Fig.  12  Nonlinear  Crack  and  Craze-Crack  Length  Behavior 


Fig.  13  Nonlinear  Behavior  of  Normalized  Craze  Length 


REMARKS  AND  SUGGESTIONS 

From  the  review  of  the  connection  between 
microstructure  and  macromechanics  as  well  as  the 
analysis  of  the  craze-crack  transition  it  is  seen  that,  in 
general,  the  complete  fracture  processes  take  time  to 
develop  and  mature.  Depending  upon  the  degree  of 
energy  absorption  by  the  microstructure  of  a  medium, 
shear  flow,  craze  or  crack  may  occur.  If  craze  occurs  first, 
then  the  transition  from  craze  to  crack  is  likely  to  be  highly 
nonlinear.  Results  are  important  in  studying  the  time- 
dependent  strength  and  fracture  behavior  of  polymeric 
and  composite  systems.  By  averaging  all  the  microscopic 
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behavior  specific  in  each  molecular  orientation,  a  fracture 
criterion  may  be  established  to  analyze  the  time  and 
deformation  dependent  breaking  strength  of  an  oriented 
polymer  solid.  Utilizing  spherical  functions  and  the 
double  Fourier  series  expansion,  the  statistical 
microscopic  behavior  in  the  vicinity  of  a  point  in  an 
oriented  medium  may  be  converted  into  several 
symmetric  tensorial  terms.  Then  the  time-to-break  can  be 
synthesized  and  incorporated  in  calculating  the 
macroscopic  behavior.  By  taking  into  account  the 
individual  values  of  the  directional  fraction  of  integrity,  the 
fracture  time  as  well  as  the  most  probable  direction  of 
fracture  initiation  within  any  volume  element  in  the 
medium  can  be  predicted.  However,  results  thus  obtained 
governs  only  the  localized  behavior  of  a  material.  For  a 
layered  composite  system,  the  interactions  among  craze- 
crack  regions  must  also  be  determined.  Perhaps  a  joint 
distribution  function  of  the  characteristic  parameters  can 
be  introduced  to  obtain  the  final  analytical  result  of  the 
strength  and  fracture  of  a  composite  system.  In  addition, 
since  the  material  systems  are  viscoelastic, 
measurements  of  temperature  fields  in  and  around  a 
craze-crack  system  are  extremely  important  in  connecting 
the  properties  of  the  microstructure  and  the  analysis  of  the 
energy  dissipation  of  the  macromechanical  benavior. 


APPENDIX 


Substituting  Equation  (42)  into  (47)  and  (48),  yields 
a 

f  $dx  *  (O2-0  j)aVc2-a2 

bJ 

+  cr,bVc2-b2  ](sin-1f-sirr1f ) 

1  _  _  aVc2-bz  -  bV  c2-a2 

+J<20r02)(a2  ♦  b2)  in 

Vc2-b2  -  V c2-af 

(2cJ,-02)ab  I 

+  2(0 |-02)a2  Inf  +  20 ,b2  Inf}.  (A1) 

c 

J ®dx  »  ^(I(o2-o  i)  aVc2-a2  ♦  o  tb^c^b^Jcos-1  f 
a 

1  ,  ,  aVcz-bz  -  bVcz-a2 

Vc2-b2  -  Vc2-a2 

+  ab°'  Vc2-“^'-a2 

+  2(OrO,)a2lnf}.  (A2) 

Differentiating  (A1)  and  (A2)  yields 


Vc2-b2  a 

+  (02-Ci)(b^=5=-  a -4a In-) 

aVc2-b5  -  bV c2-a2 


-  (02-20,)(aln 


aVc2-bs  +  bV  c2-a2 


Vc2-b2  -  V cz-a2 
'  **  \c2-b2  +  V  c2-a2  ’ 

3$ba  1,_  ,  .  <a  .  .b.c^b2 
— —  — (O  i  (sin  1  —  -  sin  1 — )  “== 


3b  * 


yfc^b T 


.  Vc2-a2  .  ..  .  b  . 

*C,,(a^P:'b  b  «  > 

aVcz-b2  -  bV  cz-az 
Vc2-b2  -  ■J  c5-a5 

-alnWrW)1*  (A4) 

~3c^  =«{[(c7r°,)VHfcz:+c,V3ferl 

/  •  .  ,b . 

(sin-1^  -  sin'1— )c 

+  (<j2'a,)  V  ‘  a'T  + 

[(02-0  ,)b2  -  o  ,a2  -  (0r20 ,)c2}} ,  (A5) 


9$ac  1.,-  _  c2-2a2  _„a  ,  .a, 

-9b"  "*  K*2'0  «><  cos  V  +  a  +  43  lnc  > 


o  .a  In 


aVc2-bz  -  bVc2-a2 
aVcz-bz  bVc2-az 


+  0,bln 


Vc2-bz  -  V c2-az 
Vc2-b2  +  V c2-a2  ’ 


3»ac  O  i  cz-2bz  a  Vc2-a2 
3b  *  *  (Ve^br 

aVc2-bz  -  bVcz-a2 
aVcz-b4  ♦  bVcz-a2 

VS**2  -  Vi21!2' » 

0  j) 4  aba 

“3?  •*{K<y2'ai)^5ipr  +a,^p!b?JCC08'1c 

-  -  .«2  -  atWc2-*2  , 

-  (o2-o,)—  -o,  Wc!;b2-}. 
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ABSTRACT 


This  paper  deals  with  a  viscoelastic  boundary  element  method  for  analyz¬ 
ing  a  polymer  quasifracture  termed  earlier  as  a  craze  in  polymers.  A  time 
dependent  boundary  stiffness  is  considered  on  the  quasifracture  envelope  sur¬ 
face.  The  viscoelastic  property  of  the  glassy  polymer  is  represented  by  a 
generalized  Kelvin  model  with  multiple  retardation  times.  According  to  the 
linear  viscoelastic  correspondence  principle,  the  associated  elasticity  solution 
can  be  solved  by  applying  the  general  integral  boundary  element  method.  Then 
the  viscoelastic  solution  in  the  time  domain  can  be  obtained  by  applying  collo¬ 
cation  Laplace  inversion  transformation.  Using  these  methods,  the  quasifrac¬ 
ture  problem  composed  of  an  isolated  craze  opening  with  time  dependent  stiff¬ 
ness  traction  in  a  stressed  rectangular  plate  is  analyzed.  The  displacement 
profile  and  the  stress  distribution  around  the  craze  envelope  surface  are  com¬ 


puted. 


2 


INTRODUCTION 

The  craze  or  quasifracture  behavior  of  glassy  polymers  has  been  studied 
by  many  scientists  using  theoretical  and/or  experimental  methods  recently 
[1.2,33.  Only  a  few  papers  reported  the  linear  elastic  quasifracture  behavior 
using  numerical  methods.  L.  Bevan  [H]  applied  both  the  elastic  finite-element 
method  and  boundary-element  method  with  linear  boundary  condition  for  investi¬ 
gating  the  craze  problem.  Recently,  using  the  nonlinear  finite  element  method, 
the  stress  distribution  around  the  envelope  surface  and  the  displacement  pro¬ 
file  associated  with  a  craze  has  been  reported  [5].  However,  since  the  boun¬ 
dary  element  method  has  currently  become  a  powerful  technique  for  solving 
boundary  value  problems  including  some  nonlinear  ones,  it  is  worthy  of  utiliz¬ 
ing,  as  it  has  several  advantages  over  the  finite  element  method.  The  number 
of  unknowns  in  the  calculating  system  depends  only  on  the  boundary  discretiza¬ 
tion  rather  than  upon  the  discretization  of  the  whole  volume  of  the  material 
body  as  in  the  finite  element  method.  The  singular  kernels  in  the  integral 
equations  weigh  the  unknown  quantities  near  a  singular  point  more  heavily  as 
compared  with  those  far  away  and  the  resultant  matrices  are  generally  well 
behaved.  The  physical  quantities  obtained  by  differentiation  of  the  primary 
variables  such  as  the  stresses  obtained  from  displacements  are  determined 
pointwise  inside  and  on  the  body.  Thus  there  is  less  chance  to  have  discontin¬ 
uities.  This  is  especially  important  in  problems  having  viscoelastic  deforma¬ 
tions  and,  in  particular,  viscoelastic  fracture  mechanics  problems  [6].  In 
addition,  this  method  takes  less  computing  time  and  yields,  greater  accuracy  as 
compared  with  those  problems  analyzed  using  the  finite  element  method  under 
somewhat  similar  situations.  Therefore,  in  the  case  wherein  highly  localized 


stresses  may  exist,  more  elements  can  be  introduced  so  that  any  possible  sin¬ 
gularities  will  not  be  suppressed  by  the  analysis.  As  the  quasifracture  beha¬ 
vior  of  crazing  is  important  for  studying  engineering  plastics  and  polymer  com¬ 
posites,  in  this  paper,  a  viscoelastic  boundary  element  method  for  analyzing 
the  polymer  quasifracture  and  determining  the  displacement  field  has  been  dev¬ 
eloped.  In  this  attempt  emphasis  is  placed  on  the  procedural  development  of 
the  method.  Actually  measured  displacement  field  obtained  earlier  has  been 
employed  [3]  in  the  computation. 

It  is  well  known  that  glassy  polymers  behave  viscoelastically.  Using  the 
correspondence  principle  in  linear  theory  oT  viscoelasticity,  the  quasifracture 
behavior  of  a  polymer  can  be  calculated  from  the  solution  of  an  associated 
elasticity  problem  by  means  of  numerical  method,  then  inversion  yields  the 
required  time  dependent  response.  Therefore,  in  this  paper  the  boundary  ele¬ 
ment  method  is  applied  to  solve  the  associated  elasticity  problem  in  the 
Laplace  domain.  By  applying  the  numerical  Laplace  inversion  technique  devel¬ 
oped  by  Schapery  [7]  and  Cost  [8],  the  associated  elasticity  solution  can  be 
transformed  from  the  Laplace  domain  back  into  the  time  domain.  There  are  sev¬ 
eral  reports  dealing  with  the  use  of  the  viscoelastic  boundary  element  method. 
For  simple  specific  viscoelastic  models,  Takashi  Kusama  and  Yasushi  Mitsui  [9] 
developed  an  improved  collocation  method  and  applied  the  boundary  element 
method  to  solve  a  Kelvin  viscoelastic  model.  Rizzo  and  Shippy  [10]  used  the 
direct  boundary  integral  method  to  solve  a" standard  linear  viscoelastic  model. 
Wang  and  Crouch  [11]  applied  the  displacement  discontinuity  boundary  element 
method  and  collocation  inversion  technique  to  solve  a  rock  mechanics  problem 
represented  by  a  Burgers  model.  In  this  paper  the  general  boundary  element 
method  together  with  the  collocation  inversion  technique  is  used  to  solve  an 
isolated  quasifracture  having  a  generalized  Kelvin  model  behavior  with  multiple 


retardation  tines.  In  using  such  a  method  the  prescribed  boundary  condition  is 
either  the  displacement  or  the  traction  condition.  For  a  quasifracture  problem 
the  boundary  condition  on  the  craze  envelope  surface  is  prescribed  in  a  stiff¬ 
ness  form.  By  considering  the  molecular  orientation  mechanism  [12]  of  the 
craze  fibril  domains,  the  boundary  displacement  of  a  craze  envelope  surface 
may  be  represented  by  a  convolution  integral.  Then  the  displacement  field  and 
the  stress  distribution  along  a  craze  surface  envelope  can  be  calculated  in 
several  time  steps.  It  is,  interesting  to  show  that  the  calculated  resulting 
stress  distribution  along  the  craze  envelope  surface  did  not  change  very  much 
with  respect  to  initial  zero  time  and  several  hundred  hours. 

FUNDAMENTAL  BOUNDARY  VALUE  PROBLEM 


The  governing  equations  for  the  quasifracture  boundary  value  problem  are 
the  equilibrium  equations  in  terms  of  the  stress  components  Ojj,  relations 
between  displacements  uj  and  strain  components  ey  together  with  a  set  of  con¬ 
stitutive  equations.  The  stress  and  displacement  fields  must  satisfy  the  pre¬ 
scribed  boundary  conditions  on  the  craze  envelope  surface  and  other  boundaries. 
They  are,  in  a  rectangular  coordinate  system  (0-Xi,x2,x,): 


The  constitutive  equations  can  be  written  in  integral  form  as: 


-  5  - 


1 


eU(t)  *  |  J(t-T)  dx. 

ekk(t)  *  |  B(t-x)  (  T)-  dx. 


(3) 


(H) 


where  J(t)  and  B(t)  are  respectively  the  creep  compliance  functions  in  shear 
and  isotropic  dilation  or  bulk  creep  compliance;  okk  and  ekk  are  respectively 
the  hydrostatic  stresses  and  strains  by  implying  the  summation  convention.  S, . 
and  e^  are  respectively  the  deviatoric  components  of  stress  tensor  Ojj  and 
strain  tensor  and  are  related  with  other  stress  and  strain  components  as 
follows: 


Sij 


°ij  ~  3  «ij°kk* 


(5) 


eij 


eij  "  3  5ijekk* 


(6) 


where  6_  are  delta  functions. 

In  a  linear  viscoelastic  polymer,  a  very  good  approximation  C3,13,1^J  for 
the  tensile  creep  compliance  D(t)  is  obtainable  using  a  generalized  Kelvin 
model  composed  of  a  series  of  Voigt  elements  or  simply  it  can  be  mathemati¬ 
cally  represented  in  the  following  form:  ' 


D(t)  = 


=  D«  ♦  dJ*1  -  exPC-~)l« 
n=i  L  nJ 


(7) 
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where  D0  and  Dn  are  constants  and  in  are  discrete  retardation  times.  Now  if 
one  adopts  the  notation  and  definition  that 

f(x,.x„s)  »  f(x,,x„t)e~St  dt,  (8) 

J  0 

where  7(x,,x»,s)  is  the  Laplace  transform  of  the  time  dependent  function 

f(xi,x„t)  with  s  as  the  Laplace  parameter,  then  it  can  be  shown  that  the  shear 
and  bulk  creep  compliance  functions  can  be  obtained  through  Laplace  inversion: 

J(t)  =  L-*(J(s»  *  L-‘[(1  <■  sv(s)U(s)]  ,  (9) 

B(t)  =  L-*(B(s))  -  L_l[( 1  -  2sv(s))B(s)]  ,  (10) 

where  u(s)  is  the  Laplace  transform  of  the  time  dependent  strain  ratio.  It  is 

to  be  noted  that  in  analyzing  a  problem  involving  the  time  dependent  viscoe¬ 
lasticity,  v(t)  is  time  dependent.  The  quantity  Poisson's  ratio  in  classical 
theory  of  elasticity  is  meaningless  in  viscoelastic  behavior,  thus  v(t)  is 

termed  strain  ratio.  Experimental  results  [15]  have  shown  that  the  strain 
ratio  v(t)  became  approximately  a  constant  for  long  creep  times.  With  this  in¬ 
formation,  the  viscoelastic  tensile  relaxation  modulus  function  E(t)  can  be 
shown  as:  __ 
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E(t)  -L'1 


1 


n 

s(D0  «• 

n=l 


Tns+1 


(n> 


In  the  craze  region,  oriented  molecular  domains  and  voids  are  formed  as 
shown  schematically  in  Fig.  1 .  Since  the  domains  are  composed  of  groups  of 
connected  fibrils  of  molecules,  they  bear  load  and  are  subjected  to  large 
deformations.  When  a  craze  elongates  its  displacement  field  in  the  direction 
of  stressing  also  increases.  While  a  part  of  the  contribution  of  the  displace¬ 
ment  field  is  because  of  the  creeping  of  the  fibril  domains,  the  major  contri¬ 
bution  comes  as  a  result  of  the  drawing  of  the  molecules  from  the  bulk  of  the 
polymer.  This  drawing  mechanism  coupled  with  a  simultaneous  neckdown  of  the 
fiber  domains  dominates  the  local  well  known  crazing  behavior  composed  of 
molecular  orientation  mechanism  and  formation  of  porosities.  A  local  strain 
field  e  (-1  <  e  <  «•  defined  as  X  -  1  where  X  is  the  draw  ratio)  identifies  the 
degree  of  molecular  orientation  termed  the  orientation  strain  and  has  been 
found  essentially  constant  [2,  3,  1 6 ]  throughout  the  craze  length  as  it  is  in¬ 
timately  associated  with  the  natural  draw  ratio  of  the  polymer.  During  the 
process  of  deformation,  the  stress  state  of  an  individual  fibril  domain  is  con¬ 
sidered  as  uniaxial  tension.  Under  a  uniaxial  stress  o«](xl,t),  a  corresponding 
small  strain  ej,(xi,t)  <<  e(Xi,t),  the  orientation  strain,  of  a  fibril  domain 
will  occur.  The  relationship  between  the  small  strain  and  the  tensile  stress 
of  each  fibril  domain  at  x,  is: 


03  3(*l,t)  - 


3cn(x„  t) 

Ef(x,,t  -  t)  - ^ -  dx. 


on  x,  <  c,  x,  =  0, 


02) 
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where  Ef(t)  is  the  viscoelastic  tensile  relaxation  modulus  of  the  fibril  domain 
while  the  orientation  strain  e  contributes  no  additional  stress.  In  prepara¬ 
tion  for  computation  using  boundary  element  method  the  traction  T,(xltt)  acting 
on  the  craze  envelope  at  x,  may  be  written  as  follows: 


Ti  (xltt) 


K(xj,t-i) 


3Uj(x,t  t) 
9t 


dt 


-  K(xlfO)  U,(x,.t) 


•t 

K(xltt-T)  U3(x„t)  di, 

.  o 


on  x,  £  c,  Xj  =  0,  (13) 


where  U3(x,,t)  is  the  opening  displacement  measured  from  the  horizontal  center 
line  of  symmetry  of  the  craze  corresponding  to  e33(x,,t)  at  the  boundary  of  the 
craze  envelope  at  x,  [2]  as  the  thickness  of  the  primordial  layer  is  small  as 
compared  with  U3.  By  writing 


K(x,,t) 


Ef(x„t) 
U3(x„0)  * 


(0  <  t  <  °>), 


(14) 


then  it  be'comes  the  stiffness  per  unit  area  of  a  craze  fibril  domain. 
K(xlt0)  =  K(x,,t) Jt=Q  is  the  initial  stiffness  at  x,.  Now  we  use  the  convolution 
integral  relationship  (13)  as  the  boundary  condition  on  the  craze  envelope. 
The  tensile  creep  compliance  Df(t)  of  a  fibril  domain  can  be  found.  By  refer¬ 
ring  to  the  molecular  orientation  theory  [12],  the  tensile  relaxation  modulus 
E f(e)  of  a  fibril  domain  may  be  represented  as  follows: 
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Ef(e)  -  CCe)E  (15) 

where,  as  stated  before,  e,  the  orientation  strain,  is  essentially  a  constant, 
thus,  C(c)  being  a  function  of  the  orientation  strain  is  also  a  constant  and  E 
is  the  modulus  of  elasticity  of  the  original  polymer  medium.  Therefore,  if  the 
time  dependency  is  introduced  as  given  in  the  following  equation: 

Ef(e,t)  -  C(c)E(t),  (16) 

both  the  nature  of  molecular  orientation  and  the  time  dependent  viscoelastic 
behavior  of  the  moduli  are  preserved.  In  the  Laplace  domain,  treat  C(e)  as 
constant,  then 


Ef(e,s)  -  C( e)E(s), 


07) 


and 


s2E(s)D(s)  =  1.  08) 

For  individual  fibril  domains  a  similar  relationship  may  be  written  as 

s^EfCsJDfis)  =  1,  (19) 

where  Df(s)  is  the  tensile  creep  compliance  function  of  the  fibril  domain  in 
the  Laplace  s-doroain.  Solving  for  Df(s),  one  gets 


Df(s)  =  C”'(  e)[s7E(s) 


(20) 


10 


or 


Df  Cs)  =  C"'(c)D(s). 


(21) 


Thus  after  inversion 


Df(t)  -  C“*( e)D(t). 


(22) 


At  position  x,t  let  us  write 


n 

Df(x»,t)  -  C",(x1.e)D(t)  =  C“(x„e){D0  ♦  £  Dn[1  -  exp(-t/T  n)] ).  (23) 

n=i 


where 


C"‘(x,,e)  =  Df  (x,,o)/D0 


(2M) 


is  a  spatial  parameter,  then  the  stiffness  becomes: 


K(xlft) 


Ef(x,,s)  K(x,,C)D0 

U|(X|.0)_  L  ^  s2D(s) 


(25) 


And  the  traction  T,(x,,t)  acting  on  the  craze  envelope  surface  equals: 


(26) 


T,(xlft)  »  L"1[sK(x,,3)U,(xl,s)  ]  . 


According  to  the  correspondence  principle  in  linear  viscoelasticity,  we 


can  transfer  the  boundary  value  problem  of  quasifracture  into  the  s-domain 
merely  by  replacing  the  elastic  parameters  by  their  corresponding  time  depen¬ 
dent  viscoelastic  parameters  in  the  Laplace  domain  s  as  follows: 


E  ■»  sl(s). 


K  -»  Kixj.OlDo/s^iXj.s), 

(Tj)q  -»  (Tj)0(s)  «  (ij)0/s, 

(Ui)o  ■»  (Ui)0(s)  -  (Ui)0/s, 


(27) 


where  (Tj)0  and  (Uj)Q  are  respectively  the  prescribed  constant  boundary  trac¬ 
tion  and  displacement  at  point  i.  Once  the  associated  elasticity  solution  is 
obtained,  then  the  Laplace  numerical  inversion  will  yield  the  time  dependent 
solution  of  the  problem. 


CALCULATING  PROCEDURE 


In  order  to  solve  the  associated  elasticity  problem,  the  general  integral 
boundary  element  method  may  be  applied.  The  detailed  investigation  of  these 
methods  and  others  may  be  found  in  the  literature  [17-193.  For  simplicity 
only  one  approach  is  utilized  3nd  the  basic  formulation  for  the  linear  elasti¬ 
city  problem  is  described  here.  In  the  two-dimensional  elastic  continuum  R 
with  boundary  f,  which  is  assumed  to  be  isotropic  without  body  force,  the  gov- 


erning  equation  may  be  obtained  from  pages  125  and  126  of  reference  [17]  as 
follows: 


ciu; 


|  T£kuk  df  =  | 


U£kTk  dr,  (k,  t,  *  1  or  3) 


(28) 


where  C*  «  1/2  for  point  i  when  it  becomes  a  boundary  point  on  a  smooth  boun- 
dary,  U£k  is  the  displacement  in  the  k  direction  due  to  a  unit  force  acting  in 
the  t  direction  at  point  i,  Uk  is  the  displacement  at  any  point  on  the  boun¬ 
dary  T  in  the  k  direction,  Tk  is  the  traction  at  any  point  on  the  boundary  r  in 
the  k  direciton,  and  T£k  is  the  traction  in  the  k  direction  due  to  a  unit  force 
acting  in  the  l  direction  at  point  i.  The  fundamental  solutions  for  the  two 
dimensional  isotropic  plane  strain  problem  are  easily  written  following  the 
equations  given  on  pages  126  and  lUl  of  reference  [,73.  They  are 


’tk 


1 


8itG(  1  -  v) 


(3  ~  ^v)6£ktn( -) 


3r 

ax£ 


Jr 

3x 


ir 

Z  • 


(29) 


'*« =  to  v)r  *Ir[n  " 2v)6**  *  2  It 


-  (1  -  2u)  nk 


3r  . 
3^  n*]}’ 


where  G,  v  are  elastic  shear  modulus  and  Poisson's  ratio  respectively  and  n^ 
is  the  outward  normal  to  the  boundary  and  r  is  the  distance  from  the  load 
point  to  the  point  under  consideration.  Equation  (29)  is  known  as  the  Kelvin's 
singular  solution  due  to  a  point  load  in  an  infinite  elastic  medium. 

At  first,  the  boundary  r  was  divided  into  N  elements  with  assumed  con¬ 
stant  values  of  Uk  and  Tk  in  each  clement.  By  applying  the  viscoelastic  cor- 
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respondence  principle,  the  following  equations  are  obtained  in  the  Laplace 
domain: 


JLf 

*  Tjkuk  dr 

q=1  J Tq 


-  ( 

)l  u£kTk  dr, 
q  =  l  ' fq 


(k,  £  =  1  or  3) 
(q  «  1,2 . N) 


(30) 


As  shown  above  there  are  2N  simultaneous  algebraic  equations.  When  2N 
boundary  tractions  and  boundary  displacements  are  given  another  2N  unknown 
boundary  tractions  and  boundary  displacements  can  be  obtained.  For  some  boun¬ 
dary  elements,  beginning  a,  say,  such  as  rq  (q  -  ct,....N),  on  which  the  stiff¬ 
ness  boundary  condition  was  prescribed.  Equation  (30)  will  become  as  follows: 


q-1 

L 


q=1 


JL  [ 

T£kUk  dr  ♦  (T£k  -  Kku£k)uk  dr  - 

rq  q*a  ' 


a- 1 

-  £  U£«Tk  dr,  (q  =  1,  2 . a . N)  (31) 

q=1  •'fq 


By  solving  the  above  simultaneous  algebraic  equations,  we  can  obtain  the 
values  Uk  and  Tk  successively  for  discrete  values  in  the  Laplace  domain. 
Based  upon  the  thermodynamic  principle,  Schapery  [7]  developed  a  collocation 
method  of  numerical  Laplace  transform  inversion.  This  method  shows  that  the 
components  of  stress  and  displacement  at  any  point  can  be  represented  by  a 
series  F(t)  defined  as  follows: 


n  - 


F(t)  =  C,  +  C2t  +  Ame~bmt, 

m  =  1 


where  Clt  C2,  Am  and  bm  are  constants.  Taking  the  Laplace  transform  of  Equa¬ 
tion  (32)  and  multiplying  by  the  transformation  parameter  s  gives: 


sF(s)  =  C 


r  m  . 

+  Si  ♦  S"  Arc 

‘  3  m'l  1  +  bm/Sj‘ 


When  time  t  goes  to  infinity  the  function  F(t)  should  be  finite.  Therefore, 
the  constant  C2  is  chosen  to  be  zero.  In  order  to  determine  the  constants  in 
this  equation  a  value  for  M  and  a  sequence  of  values  of  s  must  be  selected, 


s  =  s6  > 


(6=1,2, . ,M  +  1), 


Based  upon  Schapery's  suggestion,  the  relationship  between  s  and  t  is 
s  =  1/2t.  The  M  values  of  bm  are  taken  to  be  the  first  M+1  values  of  s.  Then 
Equation  (33)  can  be  written: 


V  An 


5F(se)  =>  Ct  ♦  — 


So 


(6  =  1,  2, . M,  M  +  1 ), 


which  is  a  set  of  M+1  linear  algebraic  equations  with  M+1  unknowns  and  Am 
solvable  using  standard  procedures.  The  guidelines  for  selecting  the  discrete 


values  of  s  can  be  found  in  Rizzo  and  Shippy  [10]. 

BEHAVIOR  OF  A  CRAZE 

According  to  the  aforementioned  theory  and  method,  an  illustration  is 
provided  by  calculating  the  displacement  field  in  the  neighborhood  of  a  hole  in 
a  linear  viscoelastic  infinite  plate.  The  load  applied  was  expressed  as  a  step 
function.  The  contour  of  the  circular  hole  has  been  divided  into  2*1  boundary 
elements  as  shown  in  Fig.  .2.  The  radial  displacements  calculated  by  either 
the  viscoelastic  boundary  element  method  or  an  analytical  solution  are  shown 
in  Fig.  3-  The  radius  of  the  hole  is  3  m.  The  applied  internal  pressure  is 
100  MN/m2.  For  linear  viscoelastic  behavior  the  tensile  creep  compliances  of 
the  material  was  represented  by  a  generalized  Kelvin  model  with  multiple  ret¬ 
ardation  times,  as  shown  below: 


=  0.238 

X 

10-3 

m2/MN, 

V 

=  0.33, 

=  0.071 

X 

10"* 

m2/MN, 

tl 

=  1 

hour, 

=  0.062 

X 

10“3 

m2/MN, 

*2 

=  10 

hours. 

=  0,0*15 

X 

10"3 

m2/MN, 

=  80 

hours, 

=  0.031 

X 

10~3 

m2/MN, 

U 

=  110 

hours. 

As  seen  in  Fig.  3.  the  computed  data  by  the  viscoelastic  boundary  element 
method  produced  excellent  agreement  with  the  analytical  results  obtained  by 
transforming  the  classical  elasticity  solution  of  a  circular  hole  in  an  infin¬ 
ite  plate  into  a  time  dependent  solution  in  linear  viscoelasticity  using  the 
well  known  correspondence  principle. 

Now  for  studying  the  quasifracture,  an  idealized  symmetrical  craze  in  a 
constant  stress  field  has  been  considered.  The  craze  basic  structure  was  rep- 
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resented  by  a  slit  with  fiber  domains  distributed  along  the  craze  envelope 
boundary.  The  distance  between  the  top  and  the  bottom  craze  envelopes  has 
been  referred  to  as  the  craze  opening  displacement  measured  from  the  center 
of  symmetry.  The  stress  acting  on  the  interface  of  the  craze  was  referred  to 
as  the  craze  envelope  stress.  The  total  craze  length  considered  was  2  mm 
thus  c  «  1  mm,  which  is  usually  referred  to  as  the  craze  length  measured  from 
the  center  of  a  craze.  The  width  of  plate  was  B  -  11.2  mm  and  the  length  of 
the  plate  was  L  =  1*t  mm  with  unit  thickness  throughout.  Because  of  symmetry, 
a  quarter.of  the  plate  containing  an  isolated  quasifracture  was  divided  into  58 
boundary  elements  as  shown  in  Fig.  it.  The  properties  of  the  bulk  material  was 
again  represented  by  a  generalized  Kelvin  model.  The  tensile  creep,  compli¬ 
ances  Dm,  retardation  times  Tm  were  given  as  before.  The  shape  of  the  applied 
stress  p0  was  a  unit  step  function  H(t)  modified  to  *42  N/mm2.  The  boundary 
conditions  used  on  the  plate  are: 


Tj(xJft)  =  0 


(38) 
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T»(x,,t)  =  0 


T,(x„t)  =  PoH(t) 


(0  S  x,  <  B,  x,  -  L).  (39) 


Uj(xj|t)  ~  0 

Ts(x,.t)  ■  0 


(x,  -  0,  0  £  x,  £  L).  (*10) 


Initially  by  using  the  finite  element  method  and  considering  the  molecular 
orientation  of  the  fibril  domains  in  the  quasifracture  [M],  the  initial  instan¬ 
taneous  craze  opening  displacement  U,(x„0)  and  the  craze  envelope  stress 
oc(x,,0)  have  been  calculated.  They  agreed  fairly  well  with  the  experimental 
results.  Subsequently  the  instantaneous  stiffness  K(x1(0)  was  calculated  and 
K(x,,t)  determined  from  Expression  (25).  In  applying  viscoelastic  boundary  ele¬ 
ment  method,  the  values  of  s-parameter  were  selected  ranging  from  10~\  10“*, 
10"*,  10°,  10\  to  1 02  and  time  t  was  chosen  as  1/25  as  given  earlier  [7-10]. 
Fig.  5  shows  the  opening  displacement  Uj  =  w  between  the  quasifracture  envel¬ 
ope  surfaces  versus  the  distance  from  the  center  of  craze  for  various  times 
corresponding  to  500,  50,  5  and  0  hours.  It  is  seen  that  the  quasifracture 
opening  displacement  increases  as  time  increases.  The  rate  of  increment  is 
relatively  high  from  0  to  50  hours.  Beyond  50  hours,  the  opening  displacement 
changes  slowly.  However,  it  is  interesting  to  find  out  that  the  stress  distri¬ 
bution  has  maintained  its  constant  value  as  shown  in  Fig.  6.  These  results  in¬ 
dicate  that  the  craze  quasifracture  behavior  can  be  successfully  analyzed 
using  this  viscoelastic  boundary  element  method. 
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FIGURE  CAPTIONS 

Fig.  1  Schemastic  Diagram  of  a  two-dimensional  craze 
Fig.  2  Mesh  division  for  viscoelastic  infinite  plate  with  a  hole 
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ABSTRACT 

In  this  paper,  the  viscoelastic  boundary  element  method  is  used  to 
estimate  the  opening  displacement  and  the  envelope  stress  on  the 
surface  of  an  isolated  crack-induced-craze  system.  To  predict  the 
propagation  history  of  both  the  crack  and. the  craze  in  a  polymer  sheet, 
the  material  properties  of  the  glassy  polymers  are  represented  by  a 
generalized  linear  viscoelastic  model.  In  the  calculation,  the  energy 
absorption  criterion  is  utilized  to  determine  the  initial  breaking  time 
and  the  propagation  rate.  A  sequence  of  numerical  calculations  of 
crack-induced-craze  propagation  by  means  of  the  boundary  element 
method  are  carried  out.  Results  are  compared  with  the  theoretical 
micromechanics  predictions.  Good  agreements  are  obtained.  This 
investigation  illustrates  that  the  three-step  envelope  stress  profile  is 
resonably  adequate  for  use  in  analyzing  polymer  quasifracture  problems. 
The  stress  concentration  phenomenon,  neglected  on  the  Dugdale  model,  is 
taken  into  consideration  in  the  present  work. 
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INTRODUCTION 


The  quasifracture  and  fracture  behaviors  of  a  crack- induced-craze 
system  in  glassy  polymers  have  been  investigated  extensively  by  many 
researchers  both  theoretical ly  and  experimentally  up  to  now1'8. 
Generally  speaking,  the  crack-induced-craze  model  in  micromechanics 
accentuats  the  opening  displacement,  the  envelope  stress  distribution 
around  the  interfaces  of  a  craze,  and  especially  the  propagation 
processes  of  the  craze  and  the  crack.  For  some  of  the  studies5"7  the 
opening  displacement  profile  has  been  determined  experimentally  first, 
from  which  the  envelope  stress  destribution  was  evaluated  by  some 
analytical  methods  such  as  Fourier  transform.  Whereas  others  obtain  the 
envelope  stresses  based  upon  experimental  observations  first,  then  the 
opening  displacement  profile  analytically.  Some  scientists6’9  considered 
that  the  yield  property  beyond  the  crack  tip  would  determine  the 
cracking  and  crazing  properties,  thus  Dugdale  model  was  employed. 
Although  Dugdale  model  is  thought  to  be  able  to  predict  the  overall 
effect,  it  is  questionable  whether  it  is  good  enough  in  representing  the 
true  behavior  of  cracking  and  crazing  properties  in  glassy  polymers. 
Based  on  a  number  of  experimental  observations  and  theoretical 
analyses6-7,  Hsiao  et  al  proposed  the  stress  step-distribution  model 
some  years  ago2,  which  seems  more  resonable  and  accurate  in 
representing  the  properties  of  the  region  behind  the  craze  tips. 

i 

Until  now  quite  a  lot  of  investigations  on  this  subject  matter  have 
been  reported,  among  which  most  were  done  by  experimental  or 


2 


analytical  method.  Papers  using  numerical  methods  for  predicting  the 
cracking  and  crazing  behaviors  have  also  appeared.  L.Bevan10,11  studied 
the  craze  micromechanics  by  using  linear  boundary  element  method,  in 
which  the  craze  at  crack  tip  was  modeled  by  linear  springs  with 
constant  stiffness.  However,  the  linear  elasticity  and  constant  stiffness 
are  not  accurate  enough  to  represent  the  properties  of 
crack -induced-craze  system  in  glassy  polymers  since  it  is  well  known 
that  glassy  polymers  behave  viscoelastically  rather  than  elastically,  and 
the  drawing  process  is  the  dominant  mechanism  in  polymer  crazing. 
Therefore,  some  scientists  considered  the  time  dependent  crack-craze 
propagation,  such  as  Chern  and  Hsiao2,  McCartney 12  and  Schapery13  who 
applied  the  linear  viscoelastic  model  for  studying  the  craze  or  crack 
propagation,  and  Schapery14.  also  studied  the  crack  growth  in 
nonhomogeneous  viscoelastic  media  for  opening  crack  model.  Some 
others15,16  investigated  the  nonlinear  quasifracture  properties  using 
finite  element  method,  and  the  time-dependent  behavior  of  a  craze  using 
viscoelastic  boundary  element  method,  in  these  studies,  the  polymer 
material  around  the  crack  or  craze  has  been  regarded  as  viscoelastic 
represented  by  a  generalized  Kelvin  model.  The  relationship  between  the 
tractions  and  the  displacement  of  fibril  domains  in  a  craze  was 
represented  by  a  convolution  integral.  Using  the  correspondence 
principle  in  linear  viscoelasticity  and  the  boundary  element  method.  :the 
time-dependent  opening  displacement  field  and  the  stress  distribution 
along  the  craze  surface  envelope  have  been  calculated  numerically. 


In  this  article,  the  viscoelastic  boundary  element  method  is  utilized 


to  study  the  propagation  of  a  crack- induced-craze  system.  Meanwhile, 
the  opening  displacement  profile  and  the  envelope  stresses  on  the  craze 
interface  surface  have  been  obtained  for  different  propagation  steps. 
For  comparison,  the  theoretical  analysis  using  energy  balance  method 
was  formulated,  with  a  three-step  stress  distribution  for  calculating  the 
propagation  of  the  crack- induced-craze  system.  And  the  opening 
displacement  profile  of  the  crack-craze  contour  was  also  evaluated. 
Because  of  the  change  of  the  boundary  conditions  during  the  propagation 
of  both  the  craze  and  the  crack,  Salamon's17  supper  posit  ion  principle  of 
a  step-like  propagation  has  been  applied  to  this  problem,  and  a  numerical 
calculation  sequence  of  the  boundary  element  method  has  been  derived. 
Comparing  the  numerical  and  analytical  results,  good  agreement  has  been 
obtained.  It  appears  that  the  step  distribution  of  the  envelope  stress 
used  in  the  analysis  is  a  good  approximation  suitable  in  dealing  with 
glassy  polymers.  The  viscoelastic  boundary  element  method  has  the 
advantage  of  ease  in  preliminary  preparation,  economical  in  computing 
time,  and  the  required  accuracy  for  studying  the  crack-induced-craze 
system  propagation  problem  may  be  achieved  without  much  difficulty. 


THE0RFT1CAI.  CONSIDERATIONS 


Craze  is  filled  with  load  bearing  highly  oriented  fibril  domains  and 
cavitated  networks  formed  by  continuous  flowing  of  the  bulk  polymer 
during  the  crazing  process.  Based  upon  some  experimental  observations 
18.19  an(j  tfie  cra2e  model  developed  earlier2,  referring  to  a  central  fixed 
(x i .  x3)  coordinate  system,  an  idealized  symmetric  crack-induced  craze 
system  in  a  constant  simple  stress  field  <?0  is  shown  in  Fig.l.  Fig. 2 
shows  c(t)  as  the  half  length  of  the  crack -craze  system  and  a(t),  the 
half  length  of  the  crack  only  at  time  t.  The  stresses  acting  on  the 
interfaces  are  called  the  envelope  stresses  with  notation  ac(x,,t)  as  a 

function  of  position  and  time.  The  half  distance  w(x,,t)  between  two 
craze  or  crack  interfaces  is  known  as  the  half  opening  displacement. 

The  half  opening  displacement  of  the  crack-induced-craze  system  in 
a  viscoelastic  polymeric  sheet  can  be  obtained  by  using  the 
correspondence  principle  in  linear  viscoelasticity  theory20.  The  field 
equations  and  the  constitutive  relations  are: 


<X„X3it)  =  0  , 

(1) 

€ij(x,,x3,t)  =  -J-lUj^X^Xj.t)  ♦  Uj.iCx^Xj.t)]  , 

(2) 

^(0  =  1  J,(t-r)— -  dr  , 

(3) 

(4) 
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where  Ojj  is  the  stress  tensor,  the  strain  tensor.  Uj  ^  the 

displacement  gradients,  and  J,  and  J2  are  respectively  the  shear  and 
isotropic  creep  compliance  functions  of  the  original  bulk  medium. 


sil<t)  =  aij(t)-U-6fjalck(t), 

(5) 

=  «!*«>, 

(6) 

with 

i.j  and  k=l,3  denoting  dummy  variables.  The  boundary 

conditions 

are. 

<?3t(x,.0.t)  =  0  .  |x,|  <c(t)  . 

(7) 

tf33(x,.0,  t)  =  tfc(x,.t)  .  |  x,  |  <c(t)  . 

(8) 

f  a33(x1.x3lt)  =  o0(t) -j 

(9) 

|  cf,|(X|,x3.t)  =  0  !■  as  (x,2  ♦  x32)->oo  . 

00) 

1  o,3(x,  ,x3.t)  =  0  J 

00 

The  opening  displacement  WQ(x,.t)  is  defined  as 

w0(x,.t)  =  U3(x,.O.t)  .  x,  <c(t)  . 

02) 

To  solve  this  viscoelastic  problem,  the  Laplace  transform  of  the 
field  equations  and  the  boundary  conditions  is  applied  to  reduce  the  time 
dependency.  Then  the  solution  to  the  transformed  equations  can  be  found 
by  using  Muskhelishvili's  complex  variable  conformal  mapping  method21 
or  Westergaard's  complex  variable  stress  function  method22  for  elastic 
medium  when  c(t),  a(t)  and  ac(x,t)  remain  unchanged  with  respect  to 
time23. 

The  time  dependent  solution  of  the  original  problem  is  obtained  by 
Laplace  inversion.  This  sofution  is  valid  only  if  the  boundary  conditions 
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are  independent  of  time  as  mentioned  earlier.  These  shortcomings  can 
be  surmounted  by  using  a  superposition  method,  i.e.  a  sequence  of  loading 
and  unloading  steps17,24,  which  yields. 

w0(x,,t)  =  Cb(0)  «(X|tt)  +  Cb(t-r)  $(x,,r)  dr  ,  (13) 

'a 

where 


4>(xlft)=cJ0‘/c2-xf  -  ^  f  1n(c+/ c2-t\2  )  dnr\ 

J  ft 


Jrf  ac(7i.t) 

*  n 


InfXf+yx,2-^2  )  dit 


wf  ac<7>’ 

JXi 


t)  lrvr(  »  ( 1 4) 


and  with  an  assumed  constant  strain  ratio  v.  l~'  being  the  Laplace 
inversion: 


Ct(t)  =  L-|-(-~|2‘,2>-  , 
s2E 


(15) 


where  E  is  the  relaxation  modulus  of  the  bulk  polymer  and  E(s) 
represents  the  same  in  the  Laplace  domain  s. 

Taking  into  consideration  of  the  thickness  of  the  primordial  layer 
from  which  the  fibril  domain  structure  has  been  pulled  out.  the  actual 
opening  displacement  of  the  system  becomes25,26: 

w(x,,t)  =  Cb(0) <J»(x,rt)  ♦  fcb(t-r)$(x,,T)dr  ♦  f  ,  (16) 

1  J0  JtX,  ' 


where  tX|  is  the  time  when  The  crack-craze-system  tip  first  reaches  the 

point  x,,  and  X  is  the  draw  ratio.  The  values  of  X  are  found  to  be 
virtually  unchanged  along  the  periphery  of  the  craze8,  27-3 '.  with  only  a 
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slightly  increase  in  the  central  region  and  near  the  craze  tip.  Therefore. 
X  can  be  considered  as  a  constant  and  the  opening  displacement  reduces 
to: 

w(x,,t)  =  -(Cb(0)4»(xt,t)  +  f  Cb(t-r)  <*>(x1tr)  dr  I , 

•'o 

for  a(t)<x,<c(t) ,  (17) 

f*.  , 

w(x1ft)  =  Cb(0)  4>(x1,t)  +  Cb(t-r)$(x1fr)  dr  +  —  w(  x,  fta)  p 

^  0 

for  0<x,<a(t)  t  (18) 

where  ta  denotes  the  time  when  the  crack  tip  arrives  at  the  point  x,=a. 

The  crack  and  craze  lengths  can  be  obtained  by  considering  the 
energy  balance.  The  energy  absorbed  by  the  craze32  is  spent  to  nucleate 
fibril  domains  near  the  craze  tips,  to  pull  fibrils  out  of  the  craze 
envelope  surfaceand  to  break  the  fibril  boundless  With  the  supposition 
that  the  fibril  nucleation  rate  at  the  craze  tips  is  proportional  to  the 
systrm  growth  velocity  c(t)33,  the  energy  rate  necessary  for  craze  tip 
growth  is  rcc(t).  Similarly,  the  energy  rate  required  for  the  crack  tip 

growth  is  raa(t),  where  rc  and  Ta  are  material  constants.  Based  upon 

the  assumption  that  the  energy,  rate  required  for  drawing  fibrils  out  of 
the  craze  envelope  surface  is  proportional  to  the  new  fibril  domain 
surface  creation  rate  S(t),  and  use  the -proportionality  constant  f  we 
have  the  following  local  energy  rate  balance  equation-’ 

Hc(t)=rcc(t)  +  raa(t)  +  rss(t)  (19) 
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w here  Hc(t)  is  the  energy  absorption  rate  of  the  quadranta!  system.  Here 

the  strain  energy,  energy  dissipation  and  kinetic  energy  have  been 
neglected  since  during  steady  state  they  are  negligible  compared  with 
the  other  rate  quantities.  With  this  simplification  and  the  terms  defined 
earlier,  it  follows  that2,3,32'34 


Hc(t)  =  |  tfc<Xi.O  — dx, , 


(20) 

(21) 


Substituting  into  (19)  results  in 


c,  „  r3Vf .  aw(x„t) „  „  .  „  . 

{csc-A~^— ) - — dx,  =  rcc  +-raa  , 


where  Vf  is  the  volume  fraction  of  the  crazed  polymer  and  d  is  the 
diameter  of  the  individual  fiber  domains. 


The  viscoelastic  boundary  element  method  is  applied  to  a  polymer 
sheet  with  an  isolated  crack- induced-craze  centrally  located.  The  sheet 
is  subjected  to  a  unit  step  tension  stress  o0H(t).  where  H(t)  is  the  unit 
step  function  of  time  t.  The  material  properties  of  the  bulk  polymer 
around  the  crack-induced-craze  system  is  considered  to  be  linearly 
viscoelastic.  The  constitutive  equation  can  be  expressed  by  convolution 
integral  equations  (3)  and  (4).  For  a  linear  viscoelastic  polymer,  a  very 
convenient  expression  for  the  creep  compliance  J  (t)  is  obtainable  by 
using  a  generalized  Kelvin  model2, 35,36  composed  of  a  seriers  of  Voigt 
elements  as  described  below: 

J(t)  =  J0  +  £  Jj(l-e 'to,  (23) 

where  J0  and  Jf  are  constants  and  Xj  ,  retardation  times.  Since  some 

experimental  results32  have  shown  that  the  strain  ratio  v(t)  remains 
virtually  unchanged  for  long  creep  times,  the  viscoelastic  relaxation 
modulus  E (t)  can  be  shown  to  be  of  the  following  form: 


E(t)=LH|i4r,=L"1 


1 


i=1 


1. 


(24) 


XjS+1 


-) 


where  the  bar  indicates  Laplace  transform  and  L'1.  Laplace  inversion. 
Taking  into  consideration  of  the  fibril  structure  of  the  craze  beyond  the 
crack  tip-"-  the  two  opposite  interfaces  of  the  craze  region  are  connected 
by  the  fibrillar  structure  as  shown  in  Fig. 2,  which  is  formed  by  a 
fibrillation  process  due  to  the  advances  of  the  crack.  These  connections 
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are  capable  of  transmitting  load  and  can  sustain  large  deformations.  The 
stress  state  of  an  individual  fibril  domain  is  considered  as  an  uniaxial 
tension.  The  relationship  between  the  traction  T3(x,.t)  and  the  opening 
displacement  U3(X|,t)  of  the  crack- induced-craze  system  can  be  shown 
as  follows1 

T3(x„0  =  -[  K(*„t-r)au-^"—  dr 

=  K(x,f0)U3(x,ft)  +  [  K(xt,t-r)  U3(xt,t)  dr  ,  (25) 

Jo 

on  Xj<c  and  x3-0, 

where  K(x,,t)  is  the  stiffness  per  unit  area  on  the  craze  surface  and 
K(Xj,0),  the  initial  value  of  IC(x,,t)  at  x,  .  Using  the  molecular 
orientation  theory16,  the  stiffness  turns  out  to  be: 

K(x„t)  =  L '[  K(SX;'^>S)J°  )  .  (26) 

It  should  be  noted  that  the  drawing  process  is  the  main  mechanism  of 
craze  thickening.  Thus  K(xt,t)  here  is  not  the  stiffness  in  the  usual 
sense.  It  must  take  the  drawing  process  into  consideration. 

The  traction  T3(x,.t)  acting  along  the  craze  envelope  becomes 

T3  (x , ,  t)  =  -  L" 1  [  S  K  (x , .  s)  U3  (x , ,  s)  ]  .  (27) 

Because  of  the  symmetry  of  the  problem,  only  a  quarter  of  the 
uniform  sheet  of  width  B,  length  L  is  considered  in  the  boundary  element 
calculation.  The  boundary  conditions  around  the  quadrantal  sheet  with 
isolated  crack -induced-craze  system  are: 


1 1 


f  T3(x„t)  =  0  , 

*  T,(x„t)  =  0  , 

0<x,sa(t),  x3=0, 

(28) 

r  T3(x,,t)  =  -  [  K(x,,t-r) 

{  Jo 

T,(x„t)  =  0  , 

9U3(x„r)  ^ 

3*  dr’ 

a(t)<x1<c(t),x3=0. 

(29) 

r  U3(X|,t)  =  0  , 
W,(x„t)=0  , 

c(t)<x1<B,  x3=0, 

(30) 

r  T,(x,,t)  =  0  , 

M3(x3,t)  =  o , 

x,=B,  0<x3<L, 

(31) 

r  Tt(Xj,t)  =  0  , 

M3(x„t)  =  <y0H(t) , 

0<x,<B,  x3=L, 

(32) 

r  U,(x3ft)  =  0  , 

W3(x3ft)  =  0, 

x,=0,  0<x3<L. 

(33) 

As  the  crack-craze  system  propagates,  new  crack  and  craze  surfaces  are 
created.  The  associated  energy  release  rate  is 

0(0  =  raa(t)  ♦  rcc(t)  .  (34) 

where  rg  and  rc  represent  respectively  the  coefficients  of  fracture  work 

for  crack  and  the  crack-induced-craze,  and  a ( t)  and  c(t)  are  respectively 
the  crack  and  the  system  propagation  rates.  The  energy  absorption  rate 
for  the  crack- induced-craze  system  is  expressed  by  Eq.  (19).  The  energy 
absorption  criterion  claims  that 
6(t)  =  Hc(t)  , 

i.e.  raa(t)  ♦  fcc(t)  =  f  ac(x|ft) dx, ,  (35) 
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In  addition,  sometimes  the  constant  crack  opening  displacement  criterion 
can  be  used  as  the  propagation  criterion  of  the  crack-induced-craze 
system.  But  these  two  criteria  are  the  same  if  the  deformation 
associated  with  the  crack  tip  is  fixed  with  respect  to  time,  a  situation 
which  occurs  when  the  applied  stress  is  constant.  Thus  both  the  crack 
and  the  craze  propagate  at  the  same  velocity37.  In  this  case,  the  opening 
displacement  of  the  crack -induced-craze  system  is  of  the  form: 

U3(x,.t)  =  flx,-c(t)]  .  (36) 

which  means  that  the  shape  of  the  crack  and  that  of  the  craze  are 
conserved  during  the  propagation,  i.e. 

^=-c(t)f'[xrc(t)l.  (3?) 

If  the  envelope  stress  is  considered  as  the  following  step  functional 
distribution: 


r 

^  «2<?0 


a(t)  <x,<b(t)  , 
b(t)  <x1<c(t)  , 


(38) 


where  c(t)-b(t)  is  a  small  quantity  representing  the  stress  growth  at  the 
craze  tip,  then  in  the  analysis,  the  energy  absorption  criterion  formula 
(35)  will  be-- 

r„a(t)  +  rcc(t)  =  f  «1a0^dx1  ♦  «2a0^dx, ,  (39) 

J  a  Jb  K' 

Tg  +  rc  =a0[cx,  U3 (a,t)-  (cx,  - o<2 )U3(b,t)] ,  (40) 

where  U3(a,t)  =  U3(x, .  t)  |  x  _  g  and  U3(b,t)  =  U3(x,.  t)|Xj  =t)  ,  are  the 
opening  displacements  of  the  crack  tip  and  the  craze  tip  respectively. 
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Noting  that  U3(b.t)<<U3(a.t). 


U3(a,t) 


ra  +  rc 


Eq.(40)  becomes- 


(41) 


which  is  the  constant  crack  opening  displacement  criterion  in  linear 
fracture  mechanics.  Therefore,  in  the  theoretical  and  numerical  analyses 
of  the  propagation  of  crack-induced-craze  system,  either  the  energy 
absorption  criterion  or  the  crack  opening  displacement  criterion  may  be 
applied  to  determine  the  propagation  rates  at  different  times. 

The  viscoelastic  boundary  element  method  for  analyzing  crack- 
induced-craze  system  in  polymers  has  been  described  in  detail  in  an 
earlier  study16.  Using  the  correspondence  principle  in  linear 
viscoelasticity  theory,  a  series  of  transformed  simultaneous  algebraic 
equations  can  be  solved.  The  displacement  Uk  and  traction  Tk  on  the 

boundary  involving  the  crack-craze  surface  can  be  obtained  in  the 
Laplace  domain.  Based  upon  Schapery’s  collection  numerical  method38 
for  Laplace  inversion,  the  components  of  the  stress  and  displacement 
fields  at  any  point  can  be  represented  for  the  fixed  time  t  by  the  series: 


F(t)  =  C0  +  C,t 


M  h  ^ 
+  £  Ame'bnlt 
m  =  l 


(42) 


with  C0 ,  C,.  Am  and  bm  being  constants.  Taking  the  Laplace 

transformation  of  Eq.  (42)  and  mutiplying  by  the  transform  parameter  s 
yield: 


sF(s)  =  C0 


m=1 


(43) 
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where-  F  (s)  designates  the  Laplace  transform  of  a  time  function  F(t). 
When  time  t  goes  to  infinity  the  function  F(t)  should  remain  finite. 
Therefore  the  constant  C,  has  to  be  chosen  as  zero.  After  a  sequence  of 
sk  (k=l,2,  ...  H+1)  is  selected,  the  constants  C0  and  Am's  can  be 
calculated  by  the  viscoelastic  boundary  element  method,  and  the  opening 
displacement  U3j  and  the  envelope  stress  T3j  on  the  jth  boundary 
etement  of  crack-craze  system  surface  for  time  t  become 

U3j=D0j+2  ,  (44) 

m  =1 

T3i=G0j-2  .  (45) 

m=1 

Substituting  the  above  two  equations  into  the  energy  absorption  criterion 
(35)  for  fracture,  the  breaking  or  the  initial  propagation  time  t  of  the 
crack-craze  system  for  discreted  boundary  elements  will  be  of  the  form: 

_  ttr(r, *  rc) 


where  iT  is  the  length  of  the  boundary  element  on  the  crack  tip  and 
the  jth  element  length  on  the  crack-craze  system,  U3j  and  T3j  are 

respectively  the  displacement  and  the  envelope  stress  on  jth  boundary 
element  of  the-  crack-craze  system  surface  before  propagation.  After 
the  commencement  of  the  propagation  of  the  crack-craze  system,  both 
the  boundary  shape  and  the  boundary  conditions  will  change  as  a  function 
of  time.  Therefore,  the  numerical  solutions  U,:  and  T^are  not  valid  for 

j  J  J 


propagating  crack-  craze  system  because  the  linear  viscoelastic 
correspondence  principle  can  only  be  applied  to  the  problem  with  time 
independent  boundary  conditions.  These  restrictions  can  be  removed  by  a 
generalized  method  of  superposition  principle,  which  uses  stepwise 
development  boundary  conditions  formulized  by  Salamon,?.  This  method 
is  utilized  here  to  deal  with  the  changing  boundary  conditions.  A  time 
dependent  function  F(t)  after  n  steps  in  the  time  interval  t  <t<tn+]  may 

be  expressed  as  follows-- 

N-1 

F(t)  =  J  (Fjfn  -  Fi[rif(t-tf<.,)]>  ♦  F„[rn  ,(t-tn)]  ,  (47) 

i=1 

where  r(  is  some  critical  linear  dimension  and  F  the  solution,  i.e.  the 
opening  displacement  or  the  traction,  which  can  be  solved  by  the  linear 
viscoelastic  correspondence  principle  in  the  ith  time  step.  tn  and  tn+ , 
are  the  nth  and  the  (n*l)th  time  steps.  Similarly,  the  envelope  traction 
T3j  and  the  displacement  U3j  on  jth  element  of  the  crack-craze  system 
after  the  nth  element  propagates  can  be  wr  iten  as: 


2  off 

i  =  1 

(g-oCjt,  e-«i(t2+t3+ — +tn) 

M 

2  off 

i  =  1 

(e-0^  -1)  +tn)  + . 

M  , 

y  o(n' 
a  ,J 

”  (e"<x»tn'1  — l)  e~oc’*n 

M 

It? 

i  =  1 

e_0<i  +  D(n) 

Oi  » 

(48) 
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T<n)  =2  Gy 
1  i=1 

.  y  G'f  _|)e-fc<t**t<*--t">  * . 

i=l  J 

♦  2  tiff*  -1)  e'^*" 

i=1 

*  2  g£’  e’O'1"  .  □<"’  .  (49) 

i=1 

According  to  the  principle,  the  calculation  procedures  are  expounded  as 
follows.  The  first  step  is  to  calculate  the  coefficients  D 4 ^  and  G 4 j on 

the“jth  element  using  the  viscoelastic  boundary  element  method,  in 
which  the  length  of  crack  is  a,  and  that  of  the  crack-craze  system  c,. 
Then  substituting  D  f  j  and  G  j  j  ( i  =  0 , 1 .  . . .  M )  into  Equation  (35)  the 

transition  or  the  initial  breaking  time  can  be  obtained.  After  that,  both 
crack  and  craze  propagate  a  distance  of  one  element  length  for  the 
steady  state  propagation  case.  Correspondingly  the  boundary  conditions 
on  the  craze  surface  will  shift  forward  by  a  length  of  one  element.  The 
second  step  is  the  calculation  of  and  (i=0,l,...M)  on  the  jth 

element  using  the  same  viscoelastic  boundary  element  method.  But  at 
this  time  the  length  of  crack  is  a,  +  £  and  that  of  system  c,+  £,  where  £ 
is  the  length  of  an  element.  Substituting  these  coefficients,  together 
with  D  -  j  and  G(j  ,  into  Equations  (25)  and  (26),  the  opening  displacement 

U3j  of  the  system  surface,  the  envelope  traction  f3)  can  be  obtained. 
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The  breaking  time  t  for  the  system  to  propagate  to  the  next  element  can 
again  be  calaulated  from  the  energy  absorption  criterion  Equation  (35). 
The  same  procedure  continues  to  be  iterated  until  the  system  grows  to 
the  nth  element.  At  that  time,  the  coefficients  D s j  and  G s -  are 

evaluated.  And  the  opening  displacement  U3j  and  the  envelope  stress  T3j 
on  the  jth  element  can  be  obtained  at  time  tn. 
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RESULTS 

A  quarter  of  the  sheet  used  in  the  calculation  using  the  boundary 
element  method  has  unit  thickness,  width  B=500pm  and  length  L=600jjm. 
The  initial  lengths  of  the  system  and  the  craze  are  taken  to  be 
c(0)=98pm  and  a(0)=38pm  respectively.  The  total  number  of  the 
boundary  elements  is  155  with  the  smallest  element  of  the  length  4pm 
located  on  the  craze  surface.  The  mesh  construction  is  shown  in  Fig. 3. 
where  the  elements  around  the  craze  tip  are  drawn  in  an  enlarged  scale. 
The  surface  of  the  crack-induced-craze  system  is  divided  into  23 
boundary  elements.  Beyond  the  craze  tip  there  are  10  elements  in  40pm 
span. 

The  viscoelastic  material  properties  are  represented  by  a  generalized 
Kelvin  model  (23)  with  other  material  constants  taken  to  be  as  follows: 
n  =  4  . 

J0  =4.17*1 0  ^  m2/MN  . 

J,  =0.71  *10~4  m2/MN  , 

J2  =0.62  x10~«  m2/MN  . 

J3  =0.43  *I0~4  m2/MN  , 

J4  =0.3  I  *10^  m2/MN  . 
v  =0.3  . 
t,  =1.0  hr.  , 
r2  =  10.0  hr.  , 
r3  =80.0  hr.  , 
r4  =  1  1 0.0  hr.  , 

Tc  =0.3  J/m2  . 
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I\  =2.72  J/m2  . 

d 

cx,  =1.21  . 

o<2  =2.58  . 

The  following  quantities  occuring  in  the  calculations  of  theoretical 
method  are  also  used39-41  = 

X  =2  . 

Ts  =0.231  J/m2  , 
d  =  4.4  nm  . 

which  represent  the  properties  of  polycarbonate.  The  applied  stress  is 
a0  =  37.4  MN/m. 

The  first  kind  of  calculation  is  based  upon  the  stiffness  distribution 
shown  in  Fig. 4.  Correspondingly,  the  instantaneous  opening  displacement 
of  the  crack-induced-craze  system  is  plotted  in  Fig. 5  against  the 
distance  measured  from  the  center  of  crack.  The  data  points  indicating 
the  experimental  observations6,8,  triangles  represent  the  theoretical 
solutions,  and  the  solid  curve  is  the  result  obtained  by  the  boundary 
element  method.  Initially  the  distribution  of  the  stress  normal  to  the 
surface  is  shown  in  Fig. 6. 

When  the  applied  constant  stress  a0  is  maintained,  the  opening 
displacement  of  the  system  increases  as  a  result  of  creep  and  the 
drawing  of  the  fibril  domains.  According  to  the  energy  absorption 
criterion,  the  craze-crack  transition  time  tb=t,  can  be  determined 

numerically.  Then  the  crack-craze  system  propagates  steadily  and  the 
case  that  the  crack  and  the  system  have  the  same  velocity9  is  considered 
here.  During  the  calculation  of  the  propagation  rate  the  stiffness  on  the 


craze -surface  is  shifted  stepwise  by  one  element.  Fig. 7  shows  the  time 
dependent  normalized  lengths  of  the  crack -induced-craze  system,  where 
the  points  are  obtained  by  the  boundary  element  method  and  the  solid 
line  is  calculated  by  theoretical  method.  As  can  be  seen  from  the  figure, 
the  propagation  rate  at  the  steady  state  is  almost  constant.  After  a 
certain  period  of  time,  both  the  crack  and  the  system  propagation  rates 
increase  drastically.  The  opening  displacement  profile  of  the 
crack-craze  system  as  a  function  of  time  is  shown  in  Fig. 8.  The  opening 
displacement  profiles  at  different  times  exhibit  somewhat  similar 
shape.  Fig.9  shows  the  comparison  between  the  results  obtained  by  the 
boundary  element  method  and  the  analytical  results  at  time  t=  12.21 
hours.  A  very  good  agreement  is  obtained.  Fig.  10  shows  the  envelope 
stress  distribution  in  the  craze  region. 

The  closeness  between  the  theoretical  results  and  those  obtained 
by  the  boundary  element  method  in  Figs.  5,  7  and  9  connotes  that  the 
boundary  element  results  are  generally  in  good  agreement  with  the 
analytical  results.  The  accuracy  depends  on  the  mesh  construction  and 
the  type  of  element  used.  Constant  elements  used  in  the  boundary 
element  calculation  procedure  yield  satisfactory  results  in  this  case. 
The  use  of  higher  order  elements,  such  as  first  and  second  order 
elements,  or  much  smaller  elements  would  improve  the  accuracy.  The 
stress  distribution  on  the  craze  surface  has  almost  the  same  shape  and 
magnitude  throughout  the  propagation.  In  fact,  it  has  been  suggested 
that  the  Dugdale  model  is  not  fully  adequate  for  analysis  in  describing 
the  craze  envelope  stress.  Nevertheless,  the  analytical  formulation 
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using,  the  three-step  distribution  function  has  been  shown  to  be  a 
resonable  and  good  approximation  for  analyzing  the  isolated  crack -craze 
system.  As  can  be  seen  in  Figs.  6  and  10.  there  is  a  deep  stress 
minimum  just  behind  the  tip  of  the  craze.  This  characteristic  feature 
persisted  during  the  course  of  this  investigation.  This  is  somewhat 
similar  to  the  results  obtained  earlier  using  the  finite  element  method. 
Two  extreme  values  in  the  envelope  stress  distribution  have  occurred. 
Like  in  the  present  case  a  minimum  envelope  stress  is  located  at  some 
point  behind  the  craze  tip  where  the  stiffness  gradient  changes  sharply 
and  a  maximum  one  occurs  at  the  tip.  -Both  of  these  extremes  have  been 
obtained^  by  either  analytical  or  experimental  methods.' •2-42-43.44  it  is 
hoped  that  this  phenomenon  will  be  studied  further  to  acquire  a  better 
understanding  of  its  behavior  with  respect  to  the  crack-craze  system. 
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FIGURES 

Fig.T  Schematic  fibrillar  structure  of  a  two-dimensional  crack - 


Fi'g.2 
Fig. 3 

Fig.d 
Fig. 5 
Fig. 6 
Fig. 7 

Fig. 8 

Fig.  9 

Fig. 10 


induced  craze  system. 

A  two-dimensional  quadrantal  crack- induced-craze  system. 
Boundary  element  mesh  with  craze  tip  shown  in  an  enlarged 
scale. 

Stiffness  distribution  in  craze  region. 

Opening  displacement  profile  of  crack-craze  system. 

Initial  envelope  stress  on  surface  of  craze  region. 

Time  dependent  normalized  length  of  crack-induced-craze 
system.  . 

Opening  displacement  of  crack-craze  system  at  several  time 
steps. 

Comparison  of  opening  displacements  by  boundary  element 
method  and  theoretical  analysis. 

Envelope  stress  on  craze  surface  at  time;  12.21  hours. 
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